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Abstract With its rapid-response capability and multiwavelength complement of instruments, 
the Swift satellite has transformed our physical understanding of gamma-ray bursts. Providing 
high-quality observations of hundreds of bursts, and facilitating a wide range of follow-up obser- 
vations within seconds of each event, Swift has revealed an unforeseen richness in observed burst 
properties, shed light on the nature of short-duration bursts, and helped realize the promise of 
gamma-ray bursts as probes of the processes and environments of star formation out to the ear- 
liest cosmic epochs. These advances have opened new perspectives on the nature and properties 
of burst central engines, interactions with the burst environment from microparsec to gigaparsec 
scales, and the possibilities for non-photonic signatures. Our understanding of these extreme 
cosmic sources has thus advanced substantially; yet more than forty years after their discovery, 
gamma-ray bursts continue to present major challenges on both observational and theoretical 
fronts. 
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I Introduction 

1.1 Setting the Stage 

Gamma-ray bursts are among the most fascinating phenomena in the universe. 
They are bright flashes of radiation with spectral energy distributions peaking 
in the 7 -ray band. They have durations measured in seconds and appear to be 
capable of producing directed flows of relativistic matter with kinetic luminosi- 
ties exceeding 10 53 erg s _1 , making them the most luminous events known. All 
evidence points to a gravitational power source associated with the cataclysmic 
formation of a relativistic star or to a precursor stage whose inevitable end point 
is a stellar mass black hole. 

The field of 7 -ray burst, alias GRB, astronomy has been greatly stimulated by the 
launch of the Swift satellite (Gehrels et al. 2004) in 2004 with its rapid response 
and panchromatic suite of instruments onboard, and by the development of new- 
technology robotic telescopes on the ground. A multidisciplinary approach is now 
emerging with data combined across the electromagnetic spectrum to learn about 
the physical processes at play; “spectral chauvinism” can no longer be tolerated in 
the modern study of GRBs. Even non-photonic neutrino and gravitational wave 
instruments are becoming more sensitive and may soon be detecting signatures 
related to GRBs. 

While interesting on their own, GRBs are now rapidly becoming powerful tools to 
study detailed properties of the galaxies in which they are embedded and of the 
universe in general. Their apparent association with massive star formation and 
their brilliant luminosities make them unique probes of the high-redshift universe 



4nd galaxy evolution. Absorption spectroscopy of GRB afterglows is being used 
to study the ISM in evolving galaxies, complementary to the traditional studies 
of quasar absorption line systems. Possibly the most interesting use of GRBs in 
cosmology is as probes of the early phases of star and galaxy formation, and the 
resulting reionization of the universe at z ~ 6 — 20. GRBs are bright enough 
to be detectable, in principle, out to much larger distances than those of the 
most luminous quasars or galaxies detected at present. Thus, promptly localized 
GRBs could serve as beacons which, shining through the pregalactic gas, provide 
information about much earlier epochs in the history of the universe. 

Before the advent of Swift, the study of GRBs had evolved somewhat unsystem- 
atically, and, as a result, the field has a great deal of historical curiosities such as 
complex classifications schemes which are now becoming streamlined as the field 
matures. Objects which were once thought to be different are now found to be 
related and the style of research has shifted from piecewise studies to a more gen- 
eral statistical approach. Although leaps in understanding can still come from 
extraordinary events as we show in several examples in this review, the appli- 
cations to broader astrophysics are coming from the compilations of hundreds 
of events. The literature on this subject has therefore become quite large, and 
we apologize for referring now and then only to the most recent comprehensive 
article in a given topic. There are several recent summary articles that give 
excellent reviews in specific areas related to GRBs. These include the supernova- 
burst connection (Woosley and Bloom 2006), short GRBs (Lee and Ramirez-Ruiz 
2007, Nakar 2007a), afterglows (van Paradijs et al. 2000, Zhang 2007) and theory 
(Meszaros 2002). Our objective here is to summarize the field of GRB astron- 



omy, from the Swift era and prior to the next steps with the Fermi Gamma Ray 
Observatory (Atwood et al. 2009), interpreting past findings while looking ahead 
to future capabilities and potential breakthroughs. 

1.2 A Burst of Progress 

The first sighting of a GRB came on July 2, 1967, from the military Vela satellites 
monitoring for nuclear explosions in verification of the Nuclear Test Ban Treaty 
(Klebesadel et al. 1973). These 7 -ray flashes, fortunately, proved to be different 
from the man-made explosions that the satellites were designed to detect and a 
new field of astrophysics was born. Over the next 30 years, hundreds of GRB 
detections were made. Frustratingly, they continued to vanish too soon to get 
an accurate angular position for follow-up observations. The reason for this is 
that 7 -rays are notoriously hard to focus, so 7 -ray images are generally not very 
sharp. 

Before 1997, most of what we knew about GRBs was based on observations from 
the Burst and Transient Source Experiment (BATSE) on board the Compton 
Gamma-Ray Observatory, whose results have been summarized in Preece et al. 
(2000). BATSE, which measured about 3000 events, revealed that between two 
and three visible bursts occur somewhere in the universe on a typical day. While 
they are on, they can outshine every other source in the 7 -ray sky, including the 
sun. Although each is unique, the bursts fall into one of two rough categories. 
Bursts that last less than two seconds are short, and those that last longer - the 
majority - are long. The two categories differ spectroscopically, with short bursts 
having relatively more high-energy 7 - rays than long bursts do. 

Arguably the most important result from BATSE concerned the distribution of 



bursts. They occur isotropically - that is, evenly over the entire sky — sug- 
gesting a cosmological distribution with no dipole and quadrupole components. 
This finding cast doubt on the prevailing wisdom, which held that bursts came 
from sources within the Milky Way. Unfortunately, 7-rays alone did not provide 
enough information to settle the question for sure. The detection of radiation 
from bursts at other wavelengths would turn out to be essential. Visible light, 
for example, could reveal the galaxies in which the bursts took place, allowing 
their distances to be measured. Attempts were made to detect these burst coun- 
terparts, but they proved fruitless. 

A watershed event occurred in 1997, when the BeppoSAX satellite succeeded 
in obtaining high-resolution X-ray images (Piro et al. 1999) of the predicted 
fading afterglow of GRB 970228 - so named because it occurred on February 28, 
1997. This detection, followed by a number of others at an approximate rate of 
10 per year, led to positions accurate to about an arc minute, which allowed the 
detection and follow-up of the afterglows at optical and longer wavelengths 1 (e.g., 
van Paradijs et al. 1997). This paved the way for the measurement of redshift 
distances, the identification of candidate host galaxies, and the confirmation that 
they were at cosmological distances (Metzger et al. 1997). 

Among the first GRBs pinpointed by BeppoSAXwas GRB 970508 (Metzger et al. 
1997). Radio observations of its afterglow provided an essential clue. The glow 
varied erratically by roughly a factor of two during the first three weeks, after 
which it stabilised and then began to diminish (Frail et al. 1997). The large varia- 
tions probably had nothing to do with the burst source itself; rather they involved 
1 We note, however, that the first optical afterglow detection of GRB 970228 (van Paradijs 

et al. 1997) was based on the X-ray prompt detection by BeppoSAX. 



the propagation of the afterglow light through space. Just as the Earth’s atrno- 
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sphere causes visible starlight to twinkle, interstellar plasma causes radio waves 
to scintillate. Therefore, if GRB 970508 was scintillating at radio wavelengths 
and then stopped, its source must have grown from a mere point to a discernible 
disk. “Discernible” here means a few light-weeks across. To reach this size the 
source must have been expanding at a considerable rate - close to the speed of 
light (Waxrnan et al. 1998). 

The observational basis for a connection between GRBs and supernovae was 
prompted by the discovery that the BeppoSAX error box of GRB 980425 con- 
tained supernova SN 1998bw (Galarna et al. 1998). A number of other GRBs 
have since shown a 1998bw-like temporal component superposed on the power 
law optical lightcurve (Woosley and Bloom 2006), but they still lacked a clear 
spectroscopic detection of an underlying supernova. Detection of such signature 
came with the discovery by the HETE-2 mission of GRB 030329 (Stanek et al. 
2003, Hjorth et al. 2003). Due to its extreme brightness and slow decay, spectro- 
scopic observations were extensive. The early spectra consisted of a power-law 
decay continuum (F u oc v~ 0 ' 9 ) typical of GRB afterglows with narrow emission 
features identifiable as Ha, [OIII], H f3 and [Oil] at z = 0.1687 (Kawabata et al. 
2003, Matheson et al. 2003) making GRB 030329 the second nearest burst over- 
all at the time and the nearest classical burst 2 . A major contribution to our 
understanding of X-ray prompt emission came also from the HETE-2 mission 
(Lamb et al. 2004) which was active from 2000 to 2006. Dozens of bursts in the 

“X-ray flash” category were observed and were found to be similar in origin to 
2 The other GRBs with supernova associations have been underluminous events (Kaneko et al. 


2007 ). 



ihe classical long GRBs (Matsuoka et al. 2004). 

Swift is the current GRB discovery mission. It is a space robot designed specif- 
ically with GRBs in mind. It combines a wide-held hard X-ray burst detection 
telescope (Burst Alert Telescope - BAT; Barthelmy et al. 2005a), with narrow- 
held X-ray (X-Ray Telescope - XRT; Burrows et al. 2005a), and UV-optical (UV 
Optical Telescope - UVOT; Rorning et al. 2004) telescopes. A powerful and fast 
on-board burst detection algorithm (Fenimore et al. 2003) provides the burst co- 
ordinates to the spacecraft, which autonomously repoints the observatory so that 
X-ray and optical observations typically commence within two minutes of the 
burst trigger. The mission was designed to hnd counterparts for all burst types, 
including the previously-elusive short GRBs. Burst positions and other data are 
provided promptly to the ground for ground observers. 

The burst detection rate for Swift is about 100 GRBs per year, resulting in a 
current data set as of December 2008 of 380 bursts. Of these, there are 126 with 
redshift determination, mostly from spectrographs on large optical telescopes and 
new robotic telescopes on the ground. These now far outnumber the ~ 40 GRB 
redshifts available prior to Swift. More than 95% of the Swift bursts have X-ray 
afterglow detection and ~ 60% have optical afterglows (UVOT + ground). To 
date 33 short bursts have been localized with 8 having redshift determinations. 
The new data have enabled much more detailed studies of the burst environment, 
the host galaxy, and the intergalactic medium between galaxies. The Swift and 
follow-up observations have also transformed our view of GRB sources. For ex- 
ample, as discussed in Section III, the old concept of a sudden release of energy 
concentrated in a few brief seconds has been discarded. Indeed, even the term 
“afterglow” is now recognized as misleading - the energy radiated during both 



phases is comparable. 
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Our primary intention in this review is to describe the most important observa- 
tional discoveries in the Swift era as well as to explain some of the new under- 
standing of key mechanisms that are believed to operate in these objects to an 
astronomical audience with little prior exposure to GRBs. Four sections follow. 
Section II provides a description of our current knowledge of what constitutes a 
GRB. Section III is a summary of the observations of the prompt and afterglow 
emission and Section IV of the observations of host galaxies and progenitor clues. 
In Section V, we examine our current progress in our understanding of the basic 
physical processes at work. Section VI is a look forward at future prospect for 
GRB study. 

II What is a Gamma-Ray Burst? 

GRBs are sudden, intense flashes of 7— rays which, for a few blinding seconds, 
light up in an otherwise fairly faint 7— ray sky. Spectra extending over many 
decades in photon energy have now been measured for hundreds of GRBs. In 
Figure 1 representative spectra are plotted in the conventional coordinates v and 
uF u , the energy radiated per logarithmic (natural log) frequency interval. Some 
obvious points should be emphasized. First, we measure directly only the energy 
radiated in the direction of the earth per second per steradian per logarithmic 
frequency interval by a source. The apparent bolometric luminosity may be quite 
different from the true bolometric luminosity if the source is not emitting isotrop- 
ically. Second, there is striking evidence for a characteristic photon energy (peak 
in the vF u spectrum), which appears to be related to the overall spectral lumi- 
nosity normalization. In contrast, the spectra of many galactic and extragalactic 



i6cretion systems are often well fitted by single power-laws. A simple power- 
law contains little information, whereas a complex spectrum composed of many 
broken power-laws tell us much more, as each break frequency must be explained. 
At cosmological distances the observed GRB fluxes imply energies that can exceed 
10 53 (O/47r) erg, where Q is the solid angle of the emitting region (Figure 2, see also 
Bloom et al. 2001). This is the mass equivalent of 0.06 Mq for the isotropic case. 
Compared with the size of the sun, the seat of this activity is extraordinarily 
compact, with sizes of less than milli-light-seconds (<300 km) as indicated by 
rapid variability of the radiation flux (Bhat et al. 1992). It is unlikely that 
mass can be converted into energy with better than a few (up to ten) percent 
efficiency; therefore, the more powerful GRB sources must “process” upwards 
of 10 _1 (fl/47r)MQ through a region which is not much larger than the size of a 
neutron star (NS) or a stellar mass black hole (BH). No other entity can convert 
mass to energy with such a high efficiency, or within such a small volume. 

The observed 7-rays have a nonthermal spectrum. Moreover, they commonly 
extend to energies above 1 MeV, the pair production threshold in the rest frame. 
These facts together imply that the emitting region must be relativistically ex- 
panding (Guilbert et al. 1983, Goodman 1986, Paczynski 1986). We draw this 
conclusion for two reasons. First, if the region were indeed only a light second 
across or less, as would be implied by the observed rapid variability in the ab- 
sence of relativistic effects, the total mass of baryons in the region would need to 
be below about 10 ~ 12 Mq in order that the electrons associated with the baryons 
should not provide a large opacity (Piran and Sherni 1993, Paczynski 1990). 
Second, larger source dimensions are required in order to avoid opacity due to 
photon-photon collisions. If the emitting region is expanding relativistically, then, 



for a given observed variation timescale, the dimension R can be increased by r 2 . 
The opacity to electrons and pairs is then reduced by T , and the threshold for 
pair production, in the observer frame, goes up by ~ T from its rest frame value 
(Fenimore et al. 1993, Woods and Loeb 1995, Baring and Harding 1997, Granot 
et al. 2008). Best-guess numbers are Lorentz factors T in the range 10 2 to 10 3 
(Lithwick and Sari 2001), allowing rapidly-variable emission to occur at radii in 
the range 10 12 to 10 14 cm. 

Because the emitting region must be several powers of ten larger than the com- 
pact object that acts as a trigger, there are further physical requirements. The 
original internal energy contained in the radiation and pairs would, after expan- 
sion, be transformed into relativistic kinetic energy. A variant that has also been 
suggested is based on the possibility that a fraction of the energy is carried by 
Poynting flux (Blandford and Znajek 1977, Usov 1992). This energy cannot be 
efficiently radiated as 7-rays unless it is re-randomized (Meszaros et al. 1994, 
Narayan et al. 1992, Rees and Meszaros 1994, Paczynski and Xu 1994). Impact 
on an external medium (or an intense external radiation field; e.g., Shaviv and 
Dar 1995) would randomize half of the initial energy merely by reducing the ex- 
pansion Lorentz factor by a factor of 2. For an approximately smooth distribution 
of external matter, the bulk Lorentz factor of the fireball thereafter decreases as 
an inverse power of the time. In the presence of turbulent magnetic fields built 
up behind the shocks (Rees and Meszaros 1992), the electrons produce a syn- 
chrotron power-law radiation spectrum which softens in time, as the synchrotron 
peak corresponding to the minimum Lorentz factor and field decreases during 
the deceleration (Katz 1994, Sari et al. 1996). Thus, the GRB radiation, which 
started out concentrated in the 7-ray range during the burst, is expected to pro- 



|?essively evolve into an afterglow radiation which peaks in the X-rays, then UV, 
optical, IR and radio (Figure 3). Detailed predictions (Meszaros and Rees 1997) 
of the afterglow properties, made in advance of the observations, agreed well 
with subsequent detections at these photon energies, followed up over periods 
of up to months. The detection of diffractive scintillation in the radio afterglow 
of GRB 970508, provided the first determination of the source size and a direct 
confirmation of relativistic source expansion (Frail et al. 1997), which was further 
strengthened by the size measurement of the afterglow image of GRB 030329 by 
radio interferometry with the VLBA (Taylor et al. 2004). 

The complex time-structure of some bursts suggests that the central engine may 
remain active for up to 100 s (Ramirez-Ruiz and Fenimore 2000) or possibly 
longer (Falcone et al. 2007). However, at much later times all memory of the 
initial time-structure would be lost: essentially all that matters is how much 
energy and momentum has been injected, its distribution in angle and velocity. 
We can at present only infer the energy per solid angle, but there are reasons 
to suspect that bursts are far from isotropic. Due to relativistic beaming, an 
observer will receive most emission from those portions of a GRB blast wave that 
are within an angle ~ 1/T of the direction to the observer. The afterglow is thus 
a signature of the geometry of the ejecta - at late stages, if the outflow is beamed, 
we expect a spherically-symmetric assumption to be inadequate; the deviations 
from the predictions of such a model would then tell us about the ejection in 
directions away from our line of sight (Rhoads 1999). 

The appearance of achromatic breaks in the development of GRB afterglows 
has been interpreted as indicating that they are jet flows beamed towards us. 
Collimation factors of Dj/47r ^ 0.01 (corresponding to half opening angles of ^ 8 



degrees) have been derived from such steepening (Frail et al. 2001, Bloom et al. 
2003). If GRB sources are beamed, then this reduces the energy per burst by two 
or three orders of magnitude at the expense of increasing their overall frequency. 
As regards the central engine trigger, there remain a number of key questions. 
What are the progenitors? What is the nature of the triggering mechanism, the 
transport of the energy and the time scales involved? Does it involve a hyper- 
accreting compact object? If so, can we decide between the various alternative 
ways of forming it? The presence of some GRBs (in the short duration category) 
in old stellar populations rules out a source uniquely associated with recent star 
formation and, in particular, massive star origin for all bursts (Gehrels et al. 
2005, Bloom et al. 2006). An understanding of the nature of these sources is thus 
inextricably linked to the metabolic pathways through which gravity, spin, and 
energy can combine to form collimated, ultrarelativistic outflows. These threads 
are few and fragile, and the tapestry is as yet a poor image of the real universe. 
If we are to improve our picture-making we must make more and stronger ties to 
physical theory. But in reconstructing the creature, we must be guided by our 
eyes and their extensions. The following sections provide a detail summary of 
the observed properties of these ultra-energetic phenomena. These threads will 
be woven in Section V. 

Ill Burst and Afterglow Observations 

The most direct diagnostics of the conditions within GRBs come from the ra- 
diation we observe, which we summarize in this section. There are many ways 
to organize a discussion of the properties of the radiation emanating from GRB 
sources. We do not intend in this article to a give a detailed review of individual 



i 4 ents since there are now sufficiently many examples that we are likely to be 
led seriously astray if we test our theories against individual events. For this 
reason, we center our discussion on major trends even in cases in which the gen- 
eralizations we describe are based on data that do not yet have the statistical 
weight of a ’’complete” sample. It should be noted that there are inherent biases 
in the discovery of a GRB at a given redshift which are often difficult to quan- 
tify, such as complex trigger efficiencies and non-detections. Continued advances 
in the observations will surely yield unexpected revisions and additions in our 
understanding of the properties of GRBs. 

III. 1 Prompt High-Energy Emission 

III. 1.1 Taxonomy The manifestations of GRB activity are extremely di- 
verse. GRBs are observed throughout the electromagnetic spectrum, from GHz 
radio waves to GeV 7-rays, but until recently, they were known predominantly 
as bursts of 7- rays, largely devoid of any observable traces at any other wave- 
lengths. 7-ray properties provide only one of several criteria for classifying GRB 
sources. Part of the problem is observational because it is not possible to obtain 
full spectral coverage in all objects and it is not easy to reconcile a classification 
based on host galaxy properties with one based on the prompt 7-ray properties. 
The major impediment to serious taxonomy is more fundamental. GRBs are het- 
erogeneous objects, especially in their directly observed properties. The success 
of a classification scheme, we believe, should be measured by the extent to which 
newly recognized properties distinguish subsets defined by differences in other 
properties. By this criterion, the taxonomy of GRBs has met with only mixed 
success. As new non-7-ray selection techniques are introduced (e.g. age of stellar 



populations in host galaxies or the presence of type Ic supernova signatures), the 
class boundaries (e.g. short and long duration events) have blurred where the 
defined subclasses transcend traditional boundaries. On the other hand, many 
new properties do correlate with old ones. This is all the more remarkable in 
that the conventional diagnostics (e.g. burst duration) measure properties on 
scales several orders of magnitude larger than that which we believe to be the 
characteristic of the powerhouse. 

III. 1.2 Observed Durations and Redshifts GRB traditionally have 
been assigned to different classes based on their duration - usually defined by 
the time during which the middle 50% (I 50 ) or 90% (T90) of the counts above 
background are measured. On the basis of this criterion, there are two classes of 
GRBs - short and long - where ~2 s duration separates them. The initial hints 
for the existence of such classes (Cline and Desai 1974, Mazets et al. 1981), were 
followed by stronger evidence from ISEE-3 and Konus-WIND data (Norris et al. 
1984) and definite proof using large statistics from BATSE (Kouveliotou et al. 
1993). BATSE results also showed that short bursts have a harder spectrum 
than long bursts (Kouveliotou et al. 1993) although this is less prominent in 
observations by Konus-WIND, HETE-2, and Swift (Sakamoto et al. 2006). 

The duration and redshift distributions for Swift GRBs are shown in Figure 4. 
The blue histogram in the left panel is for observed durations while the orange 
histogram shows the durations corrected to the source frame Tgo/(l + z) for those 
bursts with redshift determinations. In the source frame, the typical long burst 
duration is ~ 20s compared to ~ 50s in the observer frame. Swift has been 
detecting a lower fraction (~10%) of short bursts than BATSE did (25%). This 
is because Swift observes in a softer energy band (15 — 150 keV) than BATSE 



keV — 2 MeV) and because Swift requires a sky image of the event for 
burst detection and the image part of the trigger algorithm is less sensitive to 
short bursts due to their lower fluences. Shown in Figure 4 ( right panel) are the 
measured redshift distributions. The blue histogram is for Swift events while the 
gray one is for pr e-Swift, bursts. As clearly seen, Swift is currently detecting GRBs 
at a higher average redshift: < z >~ 2.5 for Swift bursts while < z >~ 1.2 for 
pr e-Swift events. The reason for this difference is the higher sensitivity of Swift 
compared to BeppoSAX and HETE-2. 

GRBs have also been classified according to their spectral properties, albeit 
less successfully. In particular, bursts with lower spectral energy peaks (-F pea k) 
have been denoted X-Ray Flashes (XRFs) based on observations by BeppoSAX , 
BATSE and HETE-2 (Heise et al. 2001, Barraud et al. 2003, Kippen et al. 2003, 
Sakamoto et al. 2005). These events are closely related to common long duration 
GRBs and appear to form a continuum of all parameters between the two types 
with no striking evidence for a discerning characteristic (Granot et al. 2005b). 

III. 1.3 Observed Correlations There is a great deal of diversity in the 
7 -ray prompt lightcurves of GRBs. Both long and short bursts can have tempo- 
ral profiles ranging from smooth, single-peaked pulses to highly structure multi- 
pulses. The prompt emission can be characterized by a variety of spectral and 
temporal parameters which include duration, variability, lag, pulse rise/fall time, 
fluence, E[ so and E pea ^. A schematic diagram illustrating the most widely dis- 
cussed 7 -ray prompt correlations is shown in Figure 5 with detailed references 
given in its caption. These correlations are often based upon statistical analysis 
of quantities whose physical causes are poorly understood, but almost certainly 
depend on many variables. Interpretations of these correlations must therefore 



be done with caution. 
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The prompt GRB lightcurves can generally be dissembled into a superposition of 
individual pulses as described by Norris et al. (1996) with rise times shorter on 
average than decay times ( panel f) . The variability or spikeness of the lightcurve 
is found to be correlated with peak luminosity or total isotropic energy of the 
burst ( panel a ) . The time lag of individual peaks seen at different energy bands 
is observed to be anticorrelated with luminosity for long bursts ( panel b ). For 
short bursts, the lag is small or not measurable. The E pe a k is also found to be 
correlated with £) so for long bursts including XRFs with short bursts as clear 
outliers ( panel c ). The total isotropic energy emission is correlated with duration 
(panel d ), with short and long bursts on approximately the same correlation line, 
albeit with a wide spread. Short bursts detected by Swift have lower E- lso on 
average than long bursts. There is a group of outliers belonging to the long burst 
category which are characterized by being significantly under luminous. These 
are GRBs: 980425, 031203, 060214. 

Numerous researchers have studied ways to determine the absolute luminosity of 
a GRB using correlations such as those illustrated in Figure 5. These include 
the lag, variability and -E pea k correlations discussed above. Other interesting cor- 
relations have included E pea ^ versus E 1 (emitted energy corrected for beaming, 
Ghirlanda et al. 2004) and E pea ^ versus a duration corrected peak luminosity 
(Firmani et al. 2006). The goal is to derive a method to determine the burst 
luminosity independently of a redshift distance determination, thus attempting 
to make GRBs standard candles which could be used, in principle, to deter- 
mine the cosmic expansion history of the universe to higher redshift than with 
supernovae. Although such efforts are currently underway (Schaefer 2006, and 



Terences therein) , it is not clear at present whether any of these correlations are 
tight enough for significant progress to be made (Bloom et al. 2003). 

III. 1.4 Soft Gamma Repeaters and Short Bursts It has been noted 
(Hurley et al. 2005, Palmer et al. 2005, Nakar et al. 2006b) that the giant flare 
(GF) observed from the putative galactic magnetar source SGR1806-20 in De- 
cember 2004 (Gaensler et al. 2005) could have looked like a typical short GRB had 
it occurred much farther away, thus making the tell-tale periodic signal charac- 
teristic of the neutron star rotation in the fading emission undetectable. The two 
previously recorded GFs of this type, one each from SGR 0520—66 on 5 March 
1979 (Fenimore et al. 1996a) and SGR 1900+14 on 27 August 1998 (Hurley et al. 
1999), would have been detectable by existing instruments only out to ~ 8 Mpc, 
and it was therefore not previously thought that they could be the source of short 
GRBs. The main spike of the 27 December event would have resembled a short, 
hard GRB if it had occurred within ~ 40 Mpc, a distance scale encompassing the 
Virgo cluster (Palmer et al. 2005). However, the paucity of observed GFs in our 
own Galaxy has so far precluded observationally based determinations of either 
their luminosity function or their rate. The observed isotropic distribution of 
short BATSE GRBs on the sky and the lack of excess events from the direction 
of the Virgo cluster suggests that only a small fraction, < 5%, of these events 
can be SGR GFs within 40 Mpc (Palmer et al. 2005). 

Before Swift detected short GRBs and their associated afterglow signatures, 
searches for nearby galaxies within narrow Inter Planetary network (IPN) er- 
ror boxes revealed that only up to ~ 15% of them could be accounted for by 
SGRs capable of producing GFs (Nakar et al. 2006b). A recent, intriguing can- 
didate is short GRB 070201 which was observed by the IPN to have a location 



consistent with the arms of the nearby (0.8 Mpc) M31 galaxy (Mazets et al. 
2008). A LIGO search for gravitational waves (Abbott et al. 2008) at the time of 
the burst turned up with no signal, thus excluding a compact merger origin. If 
the GRB was really in M31, it may have been an SGR GF. While the fraction of 
SGR events among what are now classified as short GRBs may not be dominant, 
it should be detectable and can be tested with future Swift observations. 

III. 2 Afterglow Observations 

III. 2.1 X-ray observations Swift was designed to investigate the GRB 
afterglows by filling the temporal gap between observations of the prompt emis- 
sion and the afterglow. The combined power of the BAT and XRT has revealed 
that prompt X-ray emission smoothly transitions into the decaying afterglow 
(Barthelmy et al. 2005b, O’Brien et al. 2006). Three representative Swift X-ray 
lightcurves are shown in Figure 6 for both long and short GRBs. These X-ray 
light curves start as early as a hundred seconds after the GRB trigger, and cover 
up to five decades in time. The complex behavior revealed in them sigificantly 
challenges traditional afterglow theoy, and calls into question some of the basic 
underlying assumptions. 

One of most striking results is that the early X-ray afterglows show a canonical 
behavior, where the light curve broadly consists of three distinct power law seg- 
ments (Nousek et al. 2006). A bright rapid-falling (t -a where a > 3) afterglow 
immediately after the prompt emission is followed by a steep-to-shallow transi- 
tion, which is usually accompanied by a change in the spectrum power-law index. 
This is consistent with an interpretation (Nousek et al. 2006, Zhang et al. 2006) in 
which the first break occurs when the slowly decaying emission from the forward 



Sfiock becomes dominant over the steeply decaying tail emission of the prompt 
7 -rays as seen from large angles (Kumar and Panaitescu 2000). Since these two 
components arise from physically distinct regions, their spectrum would generally 
be different. The shallow phase then transitions to the classical afterglow phase 
with no clear evidence for a spectral change. In some cases a jet break is seen at 
late times. The intermediate shallow flux stage is commonly interpreted as being 
caused by the continuous energy injection into the external shock (Nousek et al. 
2006, Zhang et al. 2006) although orientation and complex jet structures have 
been also discussed as viable alternatives. Irrespective of the underlying model, 
the presence of the shallow decay phase substantially increases the energy in the 
afterglow phase, and therefore requires the prompt 7 - 1 'ay emission mechanism to 
be very efficient. 

The average times, slopes and frequencies characterizing these three distinct X- 
ray afterglow components are listed in Table 1. Most Swift X-ray light curves are 
broadly consistent with this basic temporal description, although in most cases 
we do not see all three power law segments, either because not all are present or 
because of limited temporal coverage. The large variety of behaviors exhibited 
by afterglows at different times in their evolution can be seen in Figure 7, which 
shows the temporal history for each individual afterglow as well as the evolution 
of the cumulative X-ray afterglow luminosity for a large sample of Swift events 
with known redshift. While broadly compatible with relativistic fireball models 
(Nousek et al. 2006, Zhang et al. 2006) the complex afterglow behavior that 
has been revealed poses new challenges of interpretation. The reader is refered to 
Granot (2008) for a more detailed account of the major strengths and weaknesses 
of the standard afterglow model, as well as some of the challenges that it faces 



in explaining recent data. 
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Swift has also discovered flaring behavior appearing well after the prompt phase 
in ~ 50% of X-ray afterglow (Chincarini et al. 2007, Falcone et al. 2007). An 
illustration of bursts with bright flares is shown in Figure 6 . In some extreme 
cases, the late time flares have integrated energy similar to or exceeding the initial 
burst of 7 -rays (Burrows et al. 2005b). The rapid rise and decay, multiple flares 
in the same burst, and cases of fluence comparable to the prompt emission sug- 
gest that these flares are due to the same mechanism responsible for the prompt 
emission, which is usually attributed to the activity of the central powerhouse. 
When X-ray flares are observed by XRT, it is typically the case that no flaring is 
seen in the optical band by the UVOT. A notable example is GRB 060418 (Moli- 
nari et al. 2007), whose optical-infrared afterglow spectra is not consistent with 
a simple power-law extrapolation to soft X-ray energies and clearly requires two 
distinct spectral components. Not surprisingly, the broadband energy spectra of 
GRBs are complex and have spectra of wildly disparate shapes even for members 
belonging to a single class. 

Prior to Swift there were several reports of emission and absorption line ob- 
servations in the X-ray spectra of GRB afterglows. These included BeppoSAX 
observations of GRB 970508 (Piro et al. 1999) and GRB 000214 (Antonelli et al. 
2000), ASCA observations of GRB 970828 (Yoshida et al. 1999), Chandra obser- 
vations of GRB 991216 (Piro et al. 2000) and X MM observations of GRB 011211 
(Reeves et al. 2002). None of the detections were of high statistical significance, 
but, combined, were of some credibility. Swift has not found any significant line 
features in comprehensive observations of the X-ray afterglow of more than 200 
GRBs (Romano et al. 2008), which calls into question the significance of previous 



results. 


A final, but widely discussed, observational development is the discovery in ~ 25% 
of short bursts detected by Swift / BAT of an extended emission (EE) component 
lasting for ~100 seconds (Norris and Gehrels 2008). This component was clearly 
detected in HETE-2 burst GRB 050709 (Villasenor et al. 2005) and Swift burst 
GRB 050724 (Barthelmy et al. 2005c). Archival searches have also found BATSE 
bursts with EE (Norris and Bonnell 2006, Lazzati et al. 2001, Connaughton 2002). 
The EE is typically softer than the main peak and has an intensity range from 
10 -3 to 10 _1 times that of the initial short pulse complex. It is possible that 
many of the 75% of bursts without currently detected EE have this component 
at flux levels below detectability, although there are bursts with upper limits on 
the intensity of < 10~ 4 times that of the short pulse complex. 

III. 2. 2 Optical OBSERVATIONS With increasing frequency during the Swift 
era, optical observations of GRBs have been commencing almost immediately af- 
ter - and in one noteworthy case (GRB 080319B; Racusin et al. 2008a), prior to 
- the GRB trigger itself. These early optical observations have substantially en- 
riched our appreciation for the complexity of the physical processes active during 
the prompt emission phase of GRBs, and during the afterglow that follows. The 
large diversity of optical afterglow lightcurves can be seen in Figure 7. 

Prompt Emission. Observations prior to the Swift era demonstrated already 
that bright optical flashes such as that seen from GRB 990123 (Akerlof et al. 
1999), were relatively rare (Kehoe et al. 2001). The success of the Swift mission 
has brought a vast increase in the rate of bursts accessible to rapid optical follow- 
up, and a corresponding increase in the number of events detected from early 



times, t < 100 s. The Swift UVOT itself routinely detects optical afterglows 
following the intial prompt slew of the satellite, with a 40% detection rate for 
such bursts (Rorning et al. 2008) that is only slightly lower than the 60% detection 
rate combining all observatories. 

The fastest routine responses to Swift alerts are realized by robotic ground-based 
telescopes. The first discovery yielded by these observatories in the Swift era was 
of the 7 -ray correlated component of the prompt optical emission (Vestrand et al. 
2005, Blake et al. 2005, Vestrand et al. 2006). This component is not observed in 
every burst, but the mere fact of its correlation is sufficient to establish a common 
origin with the prompt 7 -ray emission (e.g., internal shocks). When observed, 
the ratio of the correlated 7 -ray to optical flux densities has been found to be 
roughly 10 5 to one. 

In contrast to bursts with 7 -ray correlated emission, the common burst is now 
revealed to either exhibit a single power-law decay from early times (Rykoff et al. 
2005, Quirnby et al. 2006, Yost et al. 2006), or to exhibit a flat or rising (Rykoff 
et al. 2004, Rykoff et al. 2006) or rebrightening (Stanek et al. 2007) optical 
light curve before it enters the standard power-law afterglow decay. The initial 
brightness of the typical counterpart is V ~ 14 to 17 mag (Rorning et al. 2008), 
which has made observations challenging for smaller (<1 nr) robotic facilities, and 
has limited the extent of the lightcurves collected by the UVOT. 

A few new observations of bright flaring optical emission have been collected, 
which are usually interpreted as emission from the reverse shock region (Sari 
and Piran 1999, Meszaros and Rees 1999). The early optical/NIR emission 
from GRB 041219 would have rivalled that seen from GRB 990123 if not for the 
large Galactic extinction along the line of sight (Vestrand et al. 2005). Of the 



Mree distinct peaks observed by PAIRITEL, the second may represent the on- 
set of the afterglow (reverse shock) contribution (Blake et al. 2005) - if so, a 
relatively small Lorentz factor, T ~ 70, is derived by associating the flaring 
peak time with deceleration time of the relativistic blast wave. Observations of 
GRB050525A with UVOT (Blustin et al. 2006) and GRB060111B with TAROT 
(in a unique time-resolved tracking mode; Klotz et al. 2006) show the “flattening” 
light-curve familiar from GRB 021211 (Zhang et al. 2003). Intriguingly, the “high 
redshift” GRB 050904 (Cusumano et al. 2006) at z = 6.29 (Kawai et al. 2006) 
also had prompt optical emission observed (Boer et al. 2006, Haislip et al. 2006), 
with a brightness, single-pulse structure, and fast-fading behavior reminiscent of 
GRB 990123 and thus, potentially also interpreted as reverse shock emission. 
Swift detection of the “naked eye burst” GRB 080319B, which peaked at visual 
magnitude V = 5.3 has now delivered the richest dataset, by far, addressing 
the prompt optical emission and its evolution into a standard fading afterglow 
(Racusin et al. 2008a, Bloom et al. 2009, Wozniak et al. 2008). This is only par- 
tially due to the extreme brightness of the event; the fact that it occurred at an 
equatorial location, in the night sky above the Western hemisphere, is probably 
even more important. The fact that it occurred within just one hour and ten 
degrees of the preceding GRB080319A means that it is also the only event with 
strong contraints on optical precursor emission from pointed telescopes (Racusin 
et al. 2008a). The GRB080319B dataset is rich enough that its ramifications 
are still being grappled with. In an overall sense the prompt optical emission 
correlates well with the 7 - 1 'ay lightcurve; in detail, though, the individual pulses 
observed in both bands do not track precisely, probably suggesting the presence 
of at least two distinct cooling processes at play within the dissipation region 



(Racusin et al. 2008a). 
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Afterglow Emission. The observed properties of the optical afterglow are 
largely familiar from observations prior to the Swift era (van Paradijs et al. 2000, 
Meszaros 2002) and will not be reviewed in detail here. An intriguing feature 
of later-time afterglows, revealed by the rich multiband lightcurves available in 
the Swift era, has been the occasional presence of chromatic lightcurve breaks, 
where the X-rays show a clear break (steepening of the flux decay rate) while the 
optical does not. The break in the X-ray light curve is usually identified with the 
end of the shallow decay phase. The optical lightcurve follows a single power law 
decay, usually with a temporal decay index intermediate between those in the 
X-rays before and after the break. Accommodating such chromatic breaks in a 
model where the X-ray and optical emissions arise from the same emitting region 
(Panaitescu et al. 2006) requires not only a temporal evolution of the underlying 
nricrophysical conditions within the emitting region but also fine tuning their 
photon arrival times in such a way that a break in the X-ray will be produced 
but not at optical wavelengths. Alternatively, the X-ray and optical photons 
may arise from physically distinct regions, which would naturally account for their 
seemingly decoupled behavior. Observations of the naked-eye GRB 080319B have 
sharpened this debate. This is because the optical and X-ray decays in this event 
exhibit discrepant behavior over two orders of magnitude in time, from 100 s to 
a few times 10 4 s post-burst (Racusin et al. 2008a, Kumar and Panaitescu 2008). 
At late times, t ;> 1 day, it was common in the pr e-Swift, era to observe a steepen- 
ing of the optical decay to power-law indices a Q > 2, and to associate this epoch 
with the “jet break” transition of the underlying, relativistically expanding jet. 



jet breaks are intrinsically achromatic, they are expected to manifest in the 
X-ray light curves from Swift as well as in ground-based optical observations. In 
fact, over the first three years of Swift operations only a handful of convincing 
jet-break candidates were identified (Blustin et al. 2006, Stanek et al. 2007, Dai 
et al. 2007, Willingale et al. 2007, Kocevski and Butler 2008), leading to concerns 
about the viability of this picture. 

Since then, however, deep optical imaging observations have revealed evidence 
for jet breaks in several additional Swift afterglows, a significant fraction of those 
monitored to late times (Dai et al. 2008). Systematic analysis of the Swift XRT 
data has revealed strong evidence for a jet-break signature in 12% of the X-ray 
lightcurves (Racusin et al. 2008b). A consensus has developed that jet breaks for 
typical Swift bursts may be occurring at late times and faint flux levels that are 
beyond the limit of the standard ground- and space-based campaigns. Separately, 
it has been suggested that a distinct additional spectral component could hide 
the jet-break signature in the X-ray band. For example, in the case of the bright 
GRB 070125 (Chandra et al. 2008), inverse Compton scattering of the synchrotron 
optical photons has been put forward as an explanation for the missing jet-break 
feature. 

III. 2. 3 Radio observations Radio afterglow observations are unique in 
having led to both indirect and direct demonstrations of relativistic expansion, 
via scintillation (e.g. Chandra et al. 2008) and VLBI (Taylor et al. 2005) observa- 
tions. They provide access to the properties of bright afterglows on the smallest 
angular scales. The limited number of sensitive high-resolution radio facilities, 
and the sensitivity limits of those facilities, have prevented a proportional ex- 
ploitation of the greatly-increased burst rate from Swift , primarily because the 



compensating feature of this increased burst rate has been a greater median red- 
shift and a lower characteristic afterglow flux. 

Radio afterglow observations nonetheless continue to play a vital role in accu- 
rately estimating blastwave kinetic energies (Oren et al. 2004, Granot et al. 2005a, 
Kaneko et al. 2007), with radio detections contributing crucially to the demon- 
stration of the extremely large ( E ~ 10 52 erg) kinetic energy associated with the 
high-redshift GRB 050904 (Frail et al. 2006, Gou et al. 2007), and to constrain- 
ing the relativistic energy associated with the nearby GRB 060218/SN 2006aj 
(Soderberg et al. 2006b) . Radio data can also provide a crucial “third check” on 
claims of jet-break detections, as with the broadband afterglow models applied 
to GRB 050820 (Cenko et al. 2006) and GRB 070125 (Chandra et al. 2008). 

For relatively nearby GRBs which may be associated with a supernova, radio ob- 
servations have proven invaluable in providing evidence of the physical expansion 
of GRB/SNe ejecta through direct imaging. Figure 8 shows the expansion from 
~ 0.02 mas to ~ 0.35 mas of the radio image of GRB 030329 over a time span 
from ~ 15s after the GRB trigger to 10 3 s (Pihlstrom et al. 2007). 

III. 3 Interpreting Prompt and Afterglow Emission 

The isotropic-equivalent luminosity of GRB X-ray afterglows scaled to t ~ 11 
hours after the burst in the source frame, can be used as an approximate estimator 
for the energy in the afterglow shock for the following reasons (Freedman and 
Waxrnan 2001, Piran et al. 2001, Gehrels et al. 2008). First, at 11 hr the X- 
ray band is typically above the two characteristic synchrotron frequencies, so 
that the flux has very weak dependence on nricrophysical parameters and no 
dependence on the external density, both of which are associated with relatively 



M'ge uncertainties. Second, at 11 hr the Lorentz factor of the afterglow shock 
is sufficiently small (T ~ 10) so that a large fraction of the jet is visible (out to 
an angle of 0.1 rad around the line of sight) and local inhonrogeneities 

on small angular scales are averaged out. Finally, the fact that the ratio of 
Lx( 11 hr) and Ei so is fairly constant for most GRBs, suggests that both can 
serve as a reasonable measure of the isotropic-equivalent energy content of the 
ejected outflow. 

Figure 9 shows Lx (11 hr) at 5 keV rest-frame energy as a function of their 
isotropic 7 - 1 'ay energy release for a large sample of GRBs. A linear relation, 
Lx oc (11 hr) oc L 7 i j so , seems to be broadly consistent with the data, probably 
suggesting a roughly universal efficiency for converting kinetic energy into 7 — rays 
in the prompt emission for both long and short GRBs (Lee et al. 2005a, Nousek 
et al. 2006, Bloom et al. 2007, Kaneko et al. 2007, Nysewander et al. 2008). This 
“universal” efficiency is also likely to be high (i.e., the remaining kinetic energy 
is comparable to, or even smaller than, the energy dissipated and radiated in the 
prompt emission). If this is the case, the well-known efficiency problem for long 
GRBs also persists for short events. 

III. 4 The radiated energy inventory 

With the advent of Swift, the observational inventory for GRBs has become rich 
enough to allow estimates of their radiated energy content. A compilation of the 
radiated energy in both the prompt and afterglow phases is presented in Figure 10. 
To investigate the energy dissipation behavior in the X-ray and optical afterglow, 
we fitted a natural cubic spline function to the afterglow histories (shown in 
Figure 7) for each individual afterglow and estimated the cumulative emitted 



energy. The start and end times of the integration were the first and the last 
points of the actual observations. The result is a global portrait of the effects of 
the physical processes responsible for GRB evolution, operating on scales ranging 
from AU to parsec lengths. The compilation also offers a way to assess how well 
we understand the physics of GRBs, by the degree of consistency among related 
entries. 

Figure 10 shows that the isotropic X-ray emission, Ex, iso for most Swift GRBs 
spans the range of ~ 10 50 to 10 54 erg, which is comparable to that emitted 
during the prompt 7-ray phase. Not surprisingly, events that have large isotropic- 
equivalent energy in 7- rays have larger Ex, iso; indicating a reasonably narrow 
spread in the efficiency of converting the afterglow kinetic energy into radiation. 
As can be seen also in Figure 10, the isotropic equivalent energy that is radiated 
at optical wavelengths is ~ 2 orders of magnitude smaller than that in X-rays 
and 7-rays. This is predominantly due to the fact that uF v typically peaks closer 
to the X-rays than to the optical, and it is very flat above its peak while it falls 
much faster toward lower energies. Finally, since these are isotropic equivalent 
energies, most of the contribution to the radiated afterglow energy, especially 
at radio wavelengths, is from significantly later times than for Ex, iso > and the 
collimation of the outflow together with relativistic beaming effects could result 
in much larger Ex, iso than E^\ so . 

Of the three nearby SN-GRB events plotted in Figure 10, only GRB 030329 falls 
within the ~ 10 50 erg range, while the other two events fall between 10 48 and 
~ 10 49 erg (Kaneko et al. 2007). We note that there exists a selection effect based 
on the observed photon ffux: an event is more likely to be detected when it is 
closer to us than farther, for a given intrinsic luminosity. 



@ftie of the liveliest debated issues associated with GRBs concerns the total en- 
ergy released during the explosion: are GRBs standard candles? The GRB com- 
munity has vacillated between initial claims that the GRB intrinsic luminosity 
distribution was very narrow (Horack et al. 1994) to discounting all standard can- 
dle claims, to accepting a standard total GRB energy of ~ 10 51 ergs (Frail et al. 
2001), and to diversifying GRBs into normal and sub-energetic classes (Soderberg 
et al. 2008, Ranrirez-Ruiz et al. 2005a). The important new development is that 
we now have significant observational support for the existence of a sub-energetic 
population based on the different amounts of relativistic energy released during 
the initial explosion. A network of theoretical tests lends credence to this idea 
(Waxnran 2004a, Kaneko et al. 2007, Granot and Ranrirez-Ruiz 2004, Granot 
and Kumar 2006, Pian et al. 2006). The existence of a wide range of intrinsic 
energies could further challenge the use of GRBs as standard candles. 

IV Environments and Host Galaxies 

Much of what we know about GRBs has been derived not from observations of 
the prompt burst radiation itself, but from studies of their afterglows - as they 
illuminate the circumburst surroundings - and their host galaxies. In the sections 
that follow we discuss the primary insights that have been derived in this manner 
from the study of GRB environments. 

IV. 1 Cosmological setting 

In Figure 4 we present the redshift distribution of all Swift-detected gamma-ray 
bursts. Swift and other current missions observe GRBs to cosmological distances 


quite readily; indeed, the three highest-fluence, known-redshift bursts observed 



by Swift have been at z = 0.61 (GRB 050525A), z = 2.82 (GRB 050603), and 
z = 1.26 (GRB 061007) - already spanning 40% of cosmic history. Historically, 
the majority of redshifts have been collected via host galaxy spectroscopy; in the 
Swift era, however, this pattern has been reversed - except for the short bursts 
- with the great majority of redshifts now being derived from afterglow spectra. 
In addition to theoretical arguments that posit different physical origins for short 
and long bursts (e.g., Katz and Canel 1996), the absence of short burst-associated 
supernovae to deep limits (Hjorth et al. 2005a, Fox et al. 2005, Bloom and 
Prochaska 2006, Castro-Tirado et al. 2005), 100 times fainter than SN 1998bw 
in the best cases, argues for a distinct origin of the short and long bursts. In 
agreement with this picture, the redshift distributions of the two populations are 
not consistent. 

For the long bursts, which are associated with active star formation and, in par- 
ticular, the deaths of massive stars, it is interesting to explore whether their 
distribution in redshift is consistent with other measures of cosmic star forma- 
tion. The greatly-increased number of redshifts available for Swift bursts have 
motivated several such comparisons, and for the first time, estimates of cosmic 
star formation at high redshift z > 4 using the Swift redshift sample (Chary et al. 
2007, Yiiksel et al. 2008). 

The star- formation studies over 1 < z < 4 confirm, in a broad sense, that the 
GRB redshift distribution remains consistent with independent measures of star 
formation (Jakobsson et al. 2006). However, there are signs of differential evo- 
lution of the GRB rate, in the sense that the GRB rate increases more rapidly 
with increasing redshift than expected based on star-formation measures alone 
(Le and Dernier 2007, Guetta and Piran 2007, Kistler et al. 2008). This evi- 



fiance, currently present at roughly 95%-confidence, may strengthen significantly 
in coming years. If so, this would provide a sign of bias towards low-mass and 
low-metallicity host galaxies - and potentially, low-metallicity progenitors - for 
the long-duration bursts (Section IV. 2). The theoretical curves accompanying 
the GRB redshift distribution (Figure 4) which show the evolution of a comoving 
volume element in the universe, and the volume convolved with star formation 
rate, appear to indicate that the observed distribution is wider than expected. 
The short burst redshift distribution, so far drawn exclusively from host galaxy 
observations, has also been compared to star-formation metrics. In this case, 
however, the intent has been to explore “time-delayed” progenitor models that 
correlate with star formation through a parametrized (log-normal or power-law) 
delay function (Nakar et al. 2006a, Guetta and Piran 2006, Salvaterra et al. 2008). 
Consistent with the relatively large fraction of events at z 0.5 compared to long 
bursts, and with host galaxy demographics (Zheng and Ramirez-Ruiz 2007, Shin 
and Berger 2007), these studies have concluded that a long-lived (r 1 Gyr) 
progenitor is required for these models to be consistent with Swift-e ra redshift 
measurements and the distribution of short burst fluences from BATSE. This in 
turn has led to relatively high estimates of the volumetric local short burst rate, 
at least an order of magnitude greater than the local rate of long bursts (Guetta 
et al. 2005), and a correspondingly optimistic set of predictions for Enhanced 
LIGO, VIRGO, and other ground-based gravity-wave detectors. 

IV. 2 Host galaxies of long bursts 

Surveys of GRB host galaxies in the pr e-Swift era (Le Floc’h et al. 2003, Fynbo 
et al. 2003) necessarily focused on the host galaxies of long-duration bursts. These 



surveys established a standard picture for the GRB hosts as sub-L* galaxies 
(median L ~ 0.1L*) with exponential-disk light prohles (Conselice et al. 2005, 
Wainwright et al. 2007) and high specific star- formation rates (SSFR 1 Gyr” 1 ) 
(Christensen et al. 2004). A selection of GRB host galaxies, as imaged by HST, 
are shown in Figure 11. 

This picture has not substantially changed in the Swift era. To the contrary, 
the result of various ground- and space-based efforts to characterize GRB host 
galaxies (Le Floc’h et al. 2006, Chary et al. 2007, Fynbo et al. 2008, Savaglio 
et al. 2009) have confirmed this basic outline and expanded its domain of ap- 
plicability to high redshift, combining the power of the Swift burst catalog with 
Spitzer observations. At higher redshifts, z ^ 3, it seems particularly interest- 
ing to use the GRB host galaxies to explore the evolution of mass-met allicity 
relationships that are typically compiled using held galaxies at low redshift and 
high- mass galaxy samples at high-redshift (Figure 13). Since GRB host red- 
shifts are typically secured via afterglow spectroscopy, the hosts themselves are 
uniquely free of mass and luminosity selection effects. In this context, even up- 
per limits provide useful constraints on mass-met allicity correlations and their 
evolution with redshift (Chary et al. 2007, Savaglio et al. 2009, Berger 2009). 
Such studies have also served to place GRB host galaxies in the context of other 
high-redshift galaxy populations incuding the Lyman-break and Lyman-alpha 
emitting galaxies (Fynbo et al. 2008). 

Finally, it is somewhat reassuring that numerical simulations of star- forming 
galaxy populations can generate galaxy subpopulations that reproduce the basic 
properties of the GRB host galaxies in several important respects (Courty et al. 
2004, Courty et al. 2007). In particular, selecting for high specific star- formation 



$4te generates mock galaxy catalogs with similar masses, luminosities, and colors 
to GRB host galaxies. 

IV. 2.1 Metallicity Matters While GRB host galaxies are often studied 
for the insights they provide about larger astrophysical questions, including the 
history of star formation through cosmic time, they also have the potential to 
shed light on the nature of the GRB progenitors. The association of long-duration 
bursts with star formation, for example, was proposed after observation of just 
two host galaxies (Paczynski 1998), and demonstrated firmly from the properties 
of the first twenty (Bloom et al. 2002). 

Recent years have seen a surge of interest in the question of whether GRB host 
galaxies, and hence, presumably, GRB progenitors, are nretal-poor by comparison 
to the larger population of star-forming galaxies. This question has been explored 
from a variety of perspectives, and the result of these studies, still under active 
debate, may eventually help refine our picture of the massive stellar death behind 
each GRB. 

The nretallicities of GRB host galaxies, and indeed, detailed abundance profiles, 
can be measured directly from high-resolution spectroscopy of bright afterglows. 
This area has seen dramatic progress in the Swift era, with prompt arcsecond 
positions and visual magnitude measurements from Swift UVOT and ground- 
based robotic telescopes feeding rapid-response spectroscopy from large-aperture 
facilities. One chief result of these efforts have been the collection of detailed 
abundance characterizations for multiple bursts (Vreeswijk et al. 2004, Chen et al. 
2005, Fynbo et al. 2006a, Bloom et al. 2007, Prochaska et al. 2007, Chen et al. 
2007, Price et al. 2007, Chen et al. 2008). In addition, the nretallicities of GRB 
host galaxies have been measured via standard emission-line diagnostics (Stanek 



et al. 2006, Thone et al. 2008, Savaglio et al. 2009), where possible, and in some 
cases tentative conclusions have been drawn on the basis of cruder relations such 
as mass-met allicity metrics (Berger et al. 2007a) and host galaxy luminosities 
and morphologies (Fruchter et al. 2006). 

The conclusions of these studies have yet to be reconciled into a single coherent 
picture of the nature of the GRB host galaxies and their relationship to other 
low- and high-redshift galaxy populations. However, the wealth of data do serve 
to define some associated issues with reasonable clarity. First, the metallicities of 
GRB host galaxies at z 1 are significantly ( Z ~ 0.1 Z & ) sub-solar (Savaglio et al. 
2009), consistent with the sub-solar metallicities measured for GRB host galaxies 
via absorption spectroscopy at z 2 (e.g., Chen et al. 2008). These sub-solar 
metallicities are neither surprising nor unusual for galaxy populations at high 
redshift (Fynbo et al. 2008, Savaglio et al. 2009); moreover, several candidate 
higher-metallicity hosts have been identified, although not yet confirmed (Berger 
et al. 2007a, Fruchter et al. 2006). At the same time, GRB host galaxies seem to 
be readily distinguished, in luminosity and morphology, from the host galaxies of 
core-collapse supernovae at similar redshifts (Fruchter et al. 2006), and the host 
galaxies of the lowest-redshift z <; 0.2 bursts have uniformly low metallicities that 
strongly distinguish them from the bulk of the low-redshift galaxy population 
(Stanek et al. 2006, Thone et al. 2008), and indeed, from the host galaxies of 
nearby type Ibc supernovae (Modjaz et al. 2008). 

These somewhat divergent findings might be reconciled in a picture where the 
GRB progenitors prefer (or require) a low-metallicity environment, since the in- 
creasing prevalence of such environments at z 1 would allow GRB host galaxies 
to present an increasingly fair sample of the population of star-forming galaxies 



W these higher redshifts. This argument would also dovetail with observations 
that the GRB rate seems to increase with redshift faster than the cosmic star 
formation rate, as mentioned above (Le and Dernier 2007, Guetta and Piran 
2007, Kistler et al. 2008). However, the claim that GRB host galaxies represent 
a fair sample of star-forming galaxies, even at z ^ 1, remains in dispute (e.g., 
Figure 13; Berger 2009). 

A possibly significant implication of a metallicity-dependent GRB rate would be 
an offset between the true star-formation rate and that traced by GRBs. If GRBs 
in low-metallicity environments and low-mass galaxies are more luminous, then 
they are likely over-represented in GRB samples. Low-mass galaxies and galaxy 
outskirts have lower metallicity on average and thus may yield more (and/or more 
luminous) GRBs compared to high-mass galaxies (Ramirez-Ruiz et al. 2002b). 
As galaxy mass builds up through mergers, it is also possible that the highest- 
z GRBs could be systematically more luminous due to their lower-mass host 
galaxies, an intriguing hypothesis given the extreme luminosities of some of the 
highest-redshift bursts of the Swift era, including GRB 050904. 

Finally, it is worth noting that high-resolution afterglow spectra show absorption 
features imprinted on the afterglow by gas at multiple, widely divergent physical 
scales, possibly extending from d ~ 10 pc to Gpc distances (Section IV. 2. 2). 
With the GRB metallicity question ultimately referring to the nature of the 
progenitor itself, searches for definitive signatures of the progenitor’s stellar wind 
material (Moller et al. 2002, Schaefer et al. 2003, Mirabal et al. 2003, Bloom et al. 
2007, Chen et al. 2007, Fox et al. 2008, Prochaska et al. 2008) should continue 
to be pursued and refined. At the same time, in discussing the host galaxies 
as a population, more common galaxy-integrated, emission-line diagnostics may 



better serve to place GRB hosts in the proper cosmological context. 

IV. 2. 2 Subgalactic Environments The same fast-response spectra that 
have enabled characterization of elemental abundances in GRB host galaxies have 
led to a series of discoveries regarding absorbing gas structures on sub-galactic 
scales within the GRB host: a rich array of high-ionization fine-structure Fe 
transitions in GRB 051111 (Berger et al. 2005a, Penprase et al. 2006, Prochaska 
et al. 2006); discovery of time- variability of such fine structure features in three 
bursts, demonstrating excitation by UV photons from the burst flash and young 
afterglow (Della Valle et al. 2006b, Vreeswijk et al. 2007, D’Elia et al. 2008); 
discovery of high-ionization N V features, providing evidence for absorption by 
gas within d 10 pc of the burst (Fox et al. 2008, Prochaska et al. 2008); and 
most recently, the first detection of molecular gas along the line of sight to a GRB 
afterglow (Figure 12; Prochaska et al. 2009). 

The range of these discoveries reveals a surprising complication in the interpreta- 
tion of afterglow spectra. Depending on the burst, it may be necessary to account 
for absorbing structures on every scale from the immediate parsec-size circurn- 
burst environment, to surrounding or intervening molecular clouds, to host galaxy 
absorbers on kiloparsec scales, and ultimately, to low-ionization gas associated 
with cosmological structures at gigaparsec distances. Similar complications are 
a feature of quasar absorption line studies; there, however, the situation is sim- 
plified by the persistent nature and overwhelming ionizing power of the quasar 
itself. 

Indeed, the situation may be even more complicated than this. The discovery that 
the frequency ( dN / dz ) of strong Mg II absorbers along GRB lines of sight seems 
to be roughly four times the frequency along quasar lines of sight (Prochter et al. 



§§06) suggests that some or most such absorbers seen toward GRBs - despite 
their large separation in redshift from the host galaxy - may be intrinsic to the 
GRB environment (Porciani et al. 2007). 

IV. 3 Host galaxies of short bursts 

At the time of the Swift launch, the greatest mystery in GRB astronomy was 
the origin of short GRBs. A major step forward was made in summer 2005 with 
the localization and afterglow detection of three short bursts, GRB 050509B, 
050709 and 070724. These events were found to be localized in regions of low 
star formation, either in low star forming elliptical galaxies as for GRB 050509B 
(Gehrels et al. 2005, Bloom et al. 2006) and GRB 050724 (Barthelmy et al. 2005c, 
Berger et al. 2005b) or in a region of a galaxy with low star formation (Villasenor 
et al. 2005, Fox et al. 2005, Hjorth et al. 2005b). This was in stark contrast to 
long bursts which are associated with star forming regions, and implied a non 
massive core-collapse origin. 

More than three years after these first short burst localizations, the catalog of 
confidently-identified short burst host galaxies is growing to the point where 
systematic studies can be carried out (Berger 2009) - although the effects of 
uncertain burst attribution (i.e., long or short?), uncertain host identification 
(especially for bursts with only Swift BAT or XRT localizations), and unknown 
redshifts for faint candidate hosts conspire to keep any conclusions largely quali- 
tative at this time. 

Indeed, without direct afterglow spectroscopy, association of short bursts with 
candidate host galaxies and host galaxy clusters must be approached probabilis- 
tically. In cases where a well-localized (preferably sub-arcsec) afterglow falls on a 



luminous region of the candidate host, or within a high- mass or high-redshift clus- 
ter, the association can probably be considered secure; however, in other cases an 
a posteriori estimate of the probability of association must be made (Bloom and 
Prochaska 2006, Fox and Rorning 2007). Such estimates are inevitably strongly 
dependent on input assumptions. For example: What lifetime and kick velocity 
distributions might be appropriate for progenitor binary systems? What about 
other possible progenitor classes? Any assumptions must be carefully considered 
before and after they are applied. 

With these caveats, a picture of the short burst host galaxy population as a 
whole has developed (Figure 11). It consists of three classes of host, two of which 
became apparent soon after the short burst afterglow revolution of 2005: the 
low-redshift (z 0.5), high-nrass (L ~ L*), early-type host galaxies and galaxy 
clusters, on the one hand (Gehrels et al. 2005, Fox et al. 2005, Bloom et al. 2006, 
Berger et al. 2005b, Gal- Yam et al. 2008), typified by the hosts of GRB050509B, 
GRB 050724, and GRB 050813; and the low-redshift, sub-L*, late-type galaxies 
on the other, typified by the hosts of GRB 050709 (Fox et al. 2005,Covino et al. 
2006), GRB 051221 Soderberg et al. (2006a), and GRB 061006 D’Avanzo et al. 
(2009). 

The third class of short burst host galaxies consists of faint, star- forming galaxies 
at z ^ 1 (Berger et al. 2007b, Cenko et al. 2008, Berger 2009), reminiscent 
of the host galaxies of long bursts. The existence of such higher-redshift, star- 
forming short burst host galaxies was predicted (Belczynski et al. 2006) based on 
binary population synthesis models assuming the compact object merger model 
for the short bursts. These simulations yield a bimodal distribution of lifetimes 
for merging systems, with a spike of mergers at short timescales, r ^ 100 Myr, 



followed by a dominant merger population with a r _1 lifetime distribution. 

The association of short GRBs with both star-forming and early-type galaxies has 
led to analogies with type la supernovae, whose host demographics have similarly 
provided evidence for a wide distribution of delay times between formation and 
explosion. At the same time, both core-collapse supernovae and long-duration 
GRBs are observed (almost) exclusively in late-type star-forming galaxies. As 
with supernovae and long-duration bursts, a detailed census of short burst red- 
shifts, host galaxy types, and burst locations within those hosts will undoubtedly 
help to constrain progenitor models. 

IV. 4 Neither Long Nor Short 

Some interesting gamma-ray burst host galaxies are not obviously associated with 
either the long or short burst classes. In fact, the “peculiar” cases of GRB 060505 
and GRB 060614 probably provide the first examples where studies of the host 
galaxy properties have been applied to argue for a (long- or short-) nature of the 
bursts themselves (Fynbo et al. 2006b, Gal- Yam et al. 2006, Ofek et al. 2007, 
Thone et al. 2008). 

GRB 060614 was a low-redshift, long-duration burst with no detection of a coin- 
cident supernova to deep limits. It was a bright burst (fluence in 15 — 150 keV 
band of 2.2 x 10~ 5 erg cm' 2 ) and well-studied in the X-ray and optical. With a 
T90 duration of 102 s, it seems to fall squarely in the long burst category. A host 
galaxy was found (Gal- Yam et al. 2006, Fynbo et al. 2006b, Della Valle et al. 
2006a) at z = 0.125 and deep searches were made for a coincident supernova. All 
other well-observed nearby GRBs have had supernovae detected, but this one did 
not to limits >100 times fainter than previous detections. 



41 


GRB 060614 shares some characteristics with short bursts (Gehrels et al. 2006). 
The BAT light curve shows a first short, hard-spectrum episode of emission (last- 
ing 5 s) followed by an extended and somewhat softer episode (lasting ~100 s). 
The total energy content of the second episode is 5 times that of the first. Its 
light curve shape is similar in many respects to that of short bursts with ex- 
tended emission. There are, however, differences in that the short episode of this 
event is longer than the previously detected examples and the soft episode is rel- 
atively brighter. Another similarity with short bursts comes from a lag analysis 
of GRB 060614. The lag between temporal structures in the 50 — 100 keV band 
and those in the 15 — 25 keV bands for the first 5 s is 3 ± 6 ms which falls in the 
same region of the lag-luminosity plot as short bursts (Figure 5). It is difficult 
to determine unambiguously which category of burst GRB 060614 falls into. It 
is a long event by the traditional definition, but it lacks an associated SN as 
had been observed in all other nearby long GRBs. It shares some similarities 
with Swift short bursts, but has important differences such the brightness of the 
extended soft episode. If it is due to a collapsar, it is the first indication that 
some massive star collapses either fail as supernovae (Woosley 1993) or highly 
underproduce 56 Ni (Lopez-Camara et al. 2008); if it is due to a merger, then the 
bright long-lived soft episode is hard to explain within the framework of compact 
binary mergers (Di Matteo et al. 2002, Lee et al. 2004, Setiawan et al. 2004). 
Thus, this peculiar burst has challenged the usual classifications of GRBs. 

V Basic Physical Considerations 

In this section we endeavor to outline some of the physical processes that are 
believed to be most relevant to interpreting GRBs. Though the field is far from 



lining mature, sufficient progress has been made in identifying the essential physi- 
cal ingredients. A basic scheme can provide a conceptual framework for describing 
the observations even when the framework is inaccurate! The following should 
be interpreted in this spirit. 

GRB activity manifests itself over a dynamical range of ~ 13 decades in radius. 
In Figure 11, we show a schematic montage of successive decades, exhibiting 
phenomena which are believed to take place on each of these length scales. The 
phenomena are not directly observed and the associated frames represent edu- 
cated guesses of their geometrical arrangements. An anatomical summary of the 
underlying physical processes working outward from the smallest to the largest 
scales follows. 

V.l The Central Engine 

In principle, flow onto a compact object can liberate gravitational potential en- 
ergy at a rate approaching a few tenths of Me 2 , where M is the mass inflow 
rate. Even for such high efficiencies the mass requirements are rather large, with 
the more powerful GRB sources 10 53 erg s 1 having to process upwards of 
10 _2 Mq s ' 1 through a region which is not much larger than the size of a NS 
or a stellar mass BH. Radiation of the BH rest mass on a timescale r g /c, where 
r g = GM/c 2 = 1.5 x 10 5 (M/Mq) cm is the characteristic size of the collapsed 
object, would yield luminosities c 5 /G = 4 x 10 59 erg s -1 . 

When mass accretes onto a BH or NS under these conditions, the densities and 
temperatures are so large that photons are completely trapped and neutrinos, 
being copiously emitted, are the main source of cooling. The associated interac- 
tion cross section is then many orders of magnitude smaller, and, as a result, the 



allowed accretion rates and luminosities are correspondingly much higher. For 
example, using the cross section for neutrino pair production 3 , the Eddington 
limit can be rewritten as L-^dd^ = 8 x 10 53 (E I/ /50MeV) _2 (M/M o ) erg s~ 4 , with 
an associated accretion rate is Myaa.u x (efficiency) -1 , where Mv.dd,v = LEdd,v/c 
(Ramirez-Ruiz 2006b, Lee and Ramirez-Ruiz 2007). 

The blackbody temperature if a luminosity L^dd^ emerges from a sphere of ra- 
dius 4 ?’ g is 

w V1 ~ 45 (soiv)" 72 MeV - « 

The radiation temperature is expected to be < T t h = GMm p /(3kr g ) ~ 200 MeV, 
the temperature the accreted material would reach if its gravitational potential 
energy were turned entirely into thermal energy. 

Related to this, there is a fiducial density in the vicinity of the BH 

~ 10ll(M/M °r‘ Gofsv) gcm ~ 3 - < 2) 

It should be noted that the typical Thomson optical depth under these conditions 
is tt ~ n Edd v r & ~ 10 16 and so, as expected, photons are incapable of escaping 
and constitute part of the fluid. A characteristic magnetic field strength is that 
for which B 2 /8ir = UEdd^mpC 2 : 

8 “'=(^) 1 ' , ~ J * 1 »"(‘" 1 « 1,,, ( S »lv)" G - < 3) 

If a field B^dd,u were applied to a BH with J ~ J max , the electromagnetic power 

extraction would be ~ LEdd,^ (Usov 1992). 

3 Although we have considered here the specific case of neutrino pair creation, the estimates 
vary little when one considers, for example, coherent scattering of neutrinos by nuclei and/or 
free nucleons (except for the energy scaling). 

4 Similar overall fiducial numbers also hold for neutron stars, except that the simple mass 


scalings obtained here are lost. 



^.2 Accretion Flows 


As we discussed above, a BH or NS embedded in infalling matter offers a more 
efficient power source than any other conceivable progenitor. Although this ef- 
ficiency is over a hundred times larger than that traditionally associated with 
thermonuclear reactions, the required rate of mass supply for a typical GRB is 
of course colossal. Such high mass fueling rates are never reached for BHs in 
XRBs or AGN, where the luminosity remains well below the photon Eddington 
limit. They can, however, be achieved in the process of forming NS and stellar- 
mass BHs during the collapse of massive stellar cores (Houck and Chevalier 1991, 
Woosley 1993, MacFadyen and Woosley 1999, Narayan et al. 2001) and in binary 
mergers involving compact objects (Ruffert et al. 1997, Kluzniak and Lee 1998, 
Rosswog and Liebendorfer 2003, Lee et al. 2005b, Setiawan et al. 2006, Metzger 
et al. 2008, Rosswog et al. 2004). Consequently, most recent theoretical work has 
been directed towards describing the possible formation channels for these sys- 
tems (Fryer et al. 1999, Belczynski et al. 2002, Izzard et al. 2004, Podsiadlowski 
et al. 2004). This involves evaluating those which are likely to produce a vi- 
able central engine, and understanding the flow patterns near relativistic objects 
accreting matter in the hypercritical regime (Di Matteo et al. 2002, Chen and 
Beloborodov 2007, Proga and Zhang 2006, Perna et al. 2006, Shibata et al. 2007, 
Kumar et al. 2008), where photons are unable to provide cooling but neutrinos 
do so efficiently. 

The angular momentum is quite generally a crucial parameter, in many ways 
determining the geometry of the flow (Lee and Ranrirez-Ruiz 2006). The quasi- 
spherical approximation breaks down when the gas reaches a radius r c i rc ~ 
l 2 /GM, where l is the angular momentum per unit mass, and if injection oc- 



curs more or less isotropically at large radii, a familiar accretion disk will form. 
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The inner regions of disks with mass fluxes < M Edd , v are generally able to cool 
by neutrinos on time scales shorter than the inflow time. 

If m = M/M Edd , u < 1, then the bulk of the neutrino radiation comes from 
a region only a few gravitational radii in size, and the physical conditions can 
be scaled in terms of the Eddington quantities defined above. The remaining 
relevant parameter, related to the angular momentum, is Ui n fl ow /uf re efaih where 
^freefaii — (2 GM/R) 1 ^ 2 is the free fall velocity. The inward drift speed U j n fl nw would 
be of order Vf ree f a n for supersonic radial accretion. When angular momentum is 
important, this ratio depends on the mechanism for its transport through the 
disk, which is related to the effective shear viscosity. For a thin disk, the factor 
(^inflow /'Wfreefaii) is of order a(H / R ) 2 , where H is the scale height at radius R and 
a is the phenomenological viscosity parameter (Shakura and Syunyaev 1973). 
The characteristic density, at a distance r from the hole, with account of the 
effects of rotation, is 

p ~ ?h(r/r g ) _3/2 (u inflow /Uf r eefall)7 , Edd, i /, (4) 

and the maximum magnetic field, corresponding to equipartition with the bulk 
kinetic energy, would be 

Beq ~ / (r/r g ) ^ (Tinflow/ - ^ freefall) ^ B Edd l/ . (5) 

Any neutrinos emerging directly from the central core would have energies of a 
few MeV. Note that, as mentioned above, kT Edd , v is far below the virial tem- 
perature kT v i r ~ m p c 2 (r/V g ). The flow pattern when accretion occurs would be 
then determined by the value of the parameters L p / L Edd l/ , which determine the 
importance of radiation pressure and gravity, and the ratio i C ooiAdynamicai> which 



frees the temperature if a stationary flow pattern is set up. 

It is widely believed that accretion onto a compact object, be it a NS or a stellar 
mass BH, offers the best hope of understanding the “prime mover” in all types of 
GRB sources although a possible attractive exemption includes a rapidly spinning 
neutron star with a powerful magnetic field (Usov 1992, Thompson 1994, Spruit 
1999, Thompson et al. 2004, Uzdensky and MacFadyen 2007, Bucciantini et al. 
2008, Dessart et al. 2008). If there is an ordered field B, and a characteristic 
angular velocity lo, for a spinning compact source of radius r*, then the magnetic 
dipole moment is ~ Br 3 . General arguments suggest (Gold 1968, Pacini 1967) 
that the non thermal magnetic-dipole-like luminosity will be L em ~ B 2 r®a; 4 /c 3 , 
and simple scaling from these familiar results of pulsar theory require fields of 
order 10 15 G to carry away the rotational or gravitational energy (which is ~ 
10 53 erg) in a time scale of seconds (Usov 1992). 

V.3 Jet Production 

There are two ingredients necessary for the production of jets. First, there must 
be a source of material with sufficient free energy to escape the gravitational 
field of the compact object. Second, there must be a way of imparting some 
directionality to the escaping flow. Our eventual aim must be to understand 
the overall flow pattern around a central compact object, involving accretion, 
rotation, and directional outflow but we are still far from achieving this. Most 
current workers who have discussed outflow and collimation have simply invoked 
some central supply of energy and material. A self-consistent model incorporating 
outflow and inflow must explain why some fraction of the matter can acquire a 
disproportionate share of energy (i.e., a high enthalpy). 



A spinning black hole (or neutron star) constitutes an excellent gyroscope, and 
the ingredients of accretion, angular momentum, entropy production, and possi- 
bly magnetic fields are probably sufficient to ensure the production of collimated 
outflow (McKinney 2006). However, the detailed mechanism is a matter for de- 
bate (it is not even clear what is accelerated), and several distinct lines of research 
are currently being pursued. One solution is to reconvert some of this energy via 
collisions outside the disk into electron-positron pairs or photons (Mochkovitch 
et al. 1993, Goodman et al. 1987, Rosswog and Ramirez-Ruiz 2002, Rosswog 
et al. 2003, Aloy et al. 2005, Dessart et al. 2009). If this occurs in a region of 
low baryon density (e.g., along the rotation axis, away from the equatorial plane 
of the disk) a relativistic pair-dominated wind can be produced. An obvious 
requirement for this mechanism to be efficient is that the neutrinos escape (free 
streaming, or diffusing out if the density is high enough) in a time scale shorter 
than that for advection into the BH. The efficiency for conversion into pairs (scal- 
ing with the square of the neutrino density) is too low if the neutrino production 
is too gradual, so this can become a delicate balancing act. 

Jets may alternatively be produced electromagnetically. Such mechanism could, 
in principle, circumvent the above restriction in efficiency. The potential differ- 
ence across a disk threaded by an open magnetic field lines can exceed 10 22 K, 
and this is available for accelerating high-energy particles which will produce 
an electron-positron cascade and ultimately a relativistic jet that carries away 
the binding energy of the accreting gas. Blandford & Znajek (1977) extended 
this idea to black holes and showed how the spin energy of the hole could like- 
wise be extracted. A hydromagnetic description of this mechanism is more likely 
to be appropriate (Blandford and Payne 1982). The field required to produce 



> 10 51 erg s^ 1 is enormous, and may be provided by a helical dynamo op- 
erating in hot, convective nuclear matter with a millisecond period (Duncan and 
Thompson 1992). A dipole field of the order of 10 15 G appears weak compared 
to the strongest field (Price and Rosswog 2006, Thompson et al. 2005) that can 
in principle be generated by differential rotation ( 10 17 [P/1 ms] 1 G), or by 

convection (~ 10 16 G), although how this may come about in detail is not re- 
solved. Note, however, that it only takes a residual torus (or even a cold disk) of 
10 -3 Mq to confine a field of 10 15 G. 

A potential death-trap for such relativistic outflows is the amount of entrained 
baryonic mass from the surrounding medium. For instance, a Poynting flux of 
10 53 erg could not accelerate an outflow to T > 10 2 if it had to drag more than 
~ 1O -4 M0 of baryons with it. A related complication renders the production 
of relativistic jets even more challenging, because the high neutrino fluxes are 
capable of ablating baryonic material from the surface of the disk. Thus a rest 
mass flux M r] limits the bulk Lorentz factor of the wind to T, ? = L w - m d/ (M v c 2 ) . 
Assuming that the external poloidal field strength is limited by the vigor of 
the convective motions, the spin-down luminosity scales with neutrino flux as 
^wind ~ L em oc B 2 oc u 2 on oc L^ 3 , where u con is the convective velocity. The 

pr /o 

ablation rate is oc LJ (Metzger et al. 2007, Qian and Woosley 1996), which 
indicates that the limiting bulk Lorentz factor T^ of the wind decreases as L~ l . 
Thus the burst luminosity emitted by a magnetized neutrino cooled disk may be 
self-limiting. Mass loss could, however, be suppressed if the relativistic wind were 
somehow collimated into a jet. This suggests that centrifugally driven mass loss 
will be heaviest in the outer parts of the disk, and that a detectable burst may 
be emitted only within a relatively small solid angle centered on the rotation axis 



(Levinson and Eichler 2000). 
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V.4 Jet Collimation, Stability and Confinement 

As we discussed above, one of the key issues one needs to address concerning jets 
is how they are formed? A second one is, how do they retain their coherence 
and collimation as they transverse circumburst space? The second issue is no 
less pressing than the first, since jet-like flows known on Earth are notoriously 
unstable. The situation in a jet is nonetheless different from that encountered in 
the laboratory because of the super- Alfvenic and supersonic streaming velocity 
(McKinney and Blandford 2008). MHD probably provides a better description 
of the macroscopic behavior of a GRB jet than in the case of laboratory plasmas 
because the particle Larrnor radii are so much smaller than the transverse size 
of the jet. Formal stability analysis of even the simplest jet models is complex. 
Unstable MHD modes (such as pinch and kink instability) do not grow as rapidly 
as in stationary plasma, although sufficiently short-wavelength unstable modes, 
localized within the jet interior are not suppressed by the relative motion. Lon- 
gitudinal magnetic field in the core of the jet can likewise act as a backbone, 
provided that the correlation length of the magnetic field reversal along the jet 
exceeds the wavelength of the perturbation. Instabilities may be suppressed for 
perturbations of large wavelength along the jet by the inertia of the ambient 
medium. Whatever one’s view of the relative merits of fluid and electromagnetic 
models of relativistic jets, and perhaps the truth lies between the two extremes, 
it is clear that our understanding of the role of magnetic fields in jets is less 
advanced than that of other aspects of the problem. 

An understanding of the collimation and confinement of a jet can come about 



iRly through a knowledge of the properties of the medium through which it 
propagates. Information about the ambient pressure gradient propagates into 
the jet at the internal sound speed. If the jet moves through a background with 
pressure scale height ~ r, the necessary condition for the jet interior to remain in 
pressure balance with its surroundings is 9Mj < 2, where Mj is the Mach number 
and 9 ~ d/r is the opening angle. If this condition is not satisfied when a jet is 
moving into a region of higher pressure, then strong shocks are driven into the 
jet. A jet moving into a region of lower pressure becomes overpressured relative 
to its surroundings, and thereafter expands freely. 

Freedom or confinement? This is the first question for which an extrinsic or 
environmental effect comes into play (which may, in turn, strongly affect what 
we observe), and it is thought to be particularly important for massive stellar 
GRB progenitors. This is because in such stars, a stellar envelope will remain 
to impede the advance of the jet (MacFadyen and Woosley 1999, MacFadyen 
et al. 2001, Meszaros and Rees 2001, Aloy et al. 2000, Ramirez-Ruiz et al. 2002a, 
Matzner 2003, Zhang et al. 2004, Proga et al. 2003). The beam will then evac- 
uated a channel out to some location where it impinges on the stellar envelope. 
A continuous flow of relativistic fluid emanating from the nucleus supplies this 
region with mass, momentum, energy, and magnetic flux. Most of the energy 
output during that period is deposited into a cocoon or wastebasket surround- 
ing the jet, which, after expansion, would have enough kinetic energy content to 
substantially alter the structure of the relativistic outflow, if not, in fact, provide 
much of the observed explosive power (Ramirez-Ruiz et al. 2002a). 

Can relativistic jets really be formed inside stars? This is not the most propitious 
environment to create an ultra relativistic, baryons-starved jet. What needs to be 



demonstrated, is that the outflow is not “poisoned” by baryons by the time that it 
reaches the surface of the star. It appears that it is not necessary to collimate the 
jet very tightly or to achieve high bulk Lorentz factor as the flow leaves the stellar 
surface (Zhang et al. 2004, Morsony et al. 2007). As long as the emergent flow 
has a high enthalpy per baryon, it will expand and achieve its high terminal speed 
some distance from the star (Uzdensky and MacFadyen 2007, Bucciantini et al. 
2008). A strong thermal break-out signal is expected to precede the canonical 
7 -rays observed in GRBs with massive progenitors as the shock breaks through 
the stellar surface and exposes the hot shocked material (MacFadyen et al. 2001, 
Ranrirez-Ruiz et al. 2002c, Waxnran and Meszaros 2003, Waxnran et al. 2007). 
For very extended envelopes, the jet may be unable to break through the envelope. 
TeV neutrino signals produced by Fermi accelerated relativistic protons within 
the cork may provide a means of detecting such choked-off, 7 -ray dark collapses 
(Meszaros and Waxnran 2001). 

V.5 Dissipation and Cooling Effects within the Jet 

The unique feature of GRBs is that the bulk Lorentz factor T may reach values 
from hundreds to thousands (Lithwick and Sari 2001). The relativistic motion of 
the radiating particles introduces many interesting effects in GRB emissions that 
must be properly taken into account (e.g. Feninrore et al. 1996b). 

Three frames of reference are considered when discussing the emission from sys- 
tems moving with relativistic speeds: the stationary frame, which is denoted 
here by asterisks, the comoving frame, denoted by primes, and the observer 
frame. The differential distance traveled by the expanding source during dif- 
ferential time dt* is dr = /3cdt *, where fi = \/l — T -2 . Due to time dila- 



$?:>n, dr = /3Tcdt'. The relationship between comoving and observer times is 
(1 + z)Tdt'(l — (3 cos 9) = (1 + z)dt' /5 = dt , where 0 is the angle between the 
emitting element and the observer, 6 = [T(l — /3cos#)] -1 is the Doppler factor 
and z the cosmological redshift. For an on-axis observer we therefore see that 
dt = (1 + z)dr/ r 2 c, and, as a result, the blast wave can travel a large distance 
r 2 cAf during a small observing time interval. A photon detected with dimen- 
sionless energy e = hv/m e c 2 is emitted with energy 5e’ /(I + z ). 

Few would dispute the statement that the photons which bring us all our in- 
formation about the nature of GRBs are the result of triboluminescence. For 
instance, velocity differences across the jet profile provide a source of free energy 
from particle acceleration through shock waves, hydromagnetic turbulence, and 
tearing mode magnetic reconnection (Rees and Meszaros 1994). If the value of 
r at the base increases by a factor > 2 over a timescale At, then the later ejecta 
will catch up and dissipate a fraction of their energy at radius given by 

r t ~ cAtT 2 ~ 3 x 10 14 Af 0 r 2 cm, (6) 

where At = lAto s and F = 10 2 r2- Dissipation, to be most effective, must occur 
when the wind is optically thin: tt — n'axir/ r) < 1 (here n' is the comov- 
ing number density). Otherwise it will suffer adiabatic cooling before escaping 
(Goodman 1986). The photosphere (baryonic or pair-dominated) is a source of 
soft thermalized radiation, which may be observationally detectable in some GRB 
spectra and may also result in inverse Compton cooling of the nonthermal elec- 
trons accelerated in the shocks occurring outside it (Meszaros and Rees 2000, 
Ramirez-Ruiz and Lloyd-Ronning 2002, Spada et al. 2000, Kobayashi et al. 2002, 
Meszaros et al. 2002, Ramirez-Ruiz 2005, Giannios and Spruit 2007, Pe’er et al. 
2006, Ryde and Pe’er 2008, Thompson et al. 2007). 



In the presence of turbulent magnetic fields built up behind the shocks, the 
electrons can produce a synchrotron power-law radiation spectrum, whereas the 
inverse Compton scattering of these synchrotron photons extends the spectrum 
into the GeV range (Meszaros et al. 1994). To illustrate the basic idea, suppose 
that electrons, protons, and magnetic field share the available internal energy, 
then the electrons reach typical random Lorentz factors of 7 ~ m p /m e ; while 
the assumption of a Poynting flux L b implies a comoving magnetic field of order 
B ~ Lg 2 r~ l T^ 1 ~ 10 5 ^B5o r tTi3^2 1 G, where Lb = 10 50 Lb, 50 erg s _1 and 
n = 10 13 r t ,i3- For these values of 7 and B, the typical observed synchrotron 
frequency is v sy ~ 0.5L]j ^^“^(l + z ) _1 MeV, independent of the bulk Lorentz 
factor r, and in excellent agreement with the observed values of the uF v peak 
of GRB spectra (Figure 1). Yet, there are in some instances serious problems 
associated with this model (e.g., dissipation efficiencies). These difficulties have 
motivated consideration for alternative scenarios (e.g., Kumar and Narayan 2008, 
Kumar and McMahon 2008). 

A magnetic field can ensure efficient cooling even if it is not strong enough to 
be dynamically significant. If, however, the field is dynamically significant in the 
wind (Meszaros and Rees 1997), then its internal motions could lead to dissi- 
pation even in a constant velocity wind (Thompson 1994). Instabilities in this 
magnetized wind may be responsible for particle acceleration (Thompson 2006), 
and it is possible that 7— ray production occurs mainly at large distances from 
the source (Lyutikov and Blandford 2003). 

A further effect renders the task of simulating unsteady winds even more chal- 
lenging. This stems from the likelihood that any entrained matter would be a 
mixture of protons and neutrons. If a streaming velocity builds up between ions 



§Ad neutrons then interactions can lead to dissipation even in a steady jet where 
there are no shocks (Derishev et al. 1999, Beloborodov 2003). 

V.6 Jet Interaction with the External Environment 

Astrophysicists understand supernova remnants reasonably well, despite contin- 
uing uncertainty about the initiating explosion; likewise, we may hope to under- 
stand the afterglows of GRBs, despite the uncertainties about the trigger that 
we have already emphasized. 

In the simplest version of the afterglow model, the blast wave is approximated 
by a uniform thin shell. A forward shock is formed when the expanding shell 
accelerates the external medium, and a reverse shock is formed due to deceleration 
of the cold shell. The forward and reverse shocked fluids are separated by a 
contact discontinuity and have equal kinetic energy densities. 

As the blast wave expands, it sweeps up material from the surrounding medium 
to form an external shock (Meszaros and Rees 1993). Protons captured by the 
expanding blast wave from the external medium will have total energy Tm p c 2 
in the fluid frame, where m p is the proton mass. The kinetic energy swept into 
the comoving frame by an uncollimated blast wave at the forward shock per unit 
time is given by (Blandford and McKee 1976) dE' /dt' = 47n’ 2 7i ex tm p c 3 /3r(T — 1), 
where factor of T represents the increase of external medium density due to length 
contraction, the factor (T — 1) is proportional to the kinetic energy of the swept 
up particles, and the factor /3 is proportional to the rate at which the particles 
are swept. 

The external shock becomes important when the inertia of the swept-up external 
matter starts to produce an appreciable slowing down of the ejecta. The expand- 



55 


ing shell will therefore begin to decelerate when E = TM^c 2 = r 2 m p c 2 (47rrjjn ex t/3), 
giving the deceleration radius (Meszaros and Rees 1993, Rees and Meszaros 1992) 

’' d = ( faPcS^ ) ~ 3 X 10 ‘ S (f£) 1/3 cm - <7) 

where F = E/M^c 2 is the coasting Lorentz factor, M b is the baryonic mass, 
E = 10 52 i^52 erg is the apparent isotropic energy release and n ex t is the number 
density of the circumburst medium. This sets a characteristic deceleration length. 
This deceleration allows slower ejecta to catch up, replenishing and re-energizing 
the reverse shock and boosting the momentum in the blast wave. 

Most treatments employing blast-wave theory to explain the observed afterglow 
emission from GRBs assume that the radiating particles are electrons. The prob- 
lem here is that ~ m v /m e ~ 2000 of the nonthermal particle energy swept into 
the blast-wave shock is in the form of protons or ions, unless the surroundings are 
composed primarily of electron-positron pairs. For a radiatively efficient system, 
physical processes must therefore transfer a large fraction of the swept-up energy 
to the electron component (Gedalin et al. 2008). In most elementary treatments 
it is simply assumed that a fraction e e of the forward-shock power is transferred 
to the electrons. 

The strength of the magnetic field is another major uncertainty. The stan- 
dard prescription is to assume that the magnetic field energy density mb is a 
fixed fraction e# of the downstream energy density of the shocked fluid. Hence 
ub = B 2 /(8n) = 4esn ex tm p c 2 (r 2 — T) (although see, e.g., Rossi and Rees 
2003). It is also generally supposed in simple blast- wave model calculations that 
some mechanism injects electrons with a power-law distribution between elec- 
tron Lorentz factors 7min < 7 < 7max downstream of the shock front, where the 
maximum injection energy is obtained by balancing synchrotron losses and an 



iiceleration rate given in terms of the inverse of the Larrnor time scale. 

A break is formed in the electron spectrum at cooling electron Lorentz factor 
7c, which is found by balancing the synchrotron loss time scale t' sy with the adi- 


abatic expansion time t' adi ~ r/(Tc) (Sari et al. 1998). For an adiabatic blast 
wave, T oc f _3//8 , so that 7 m ; n oc f -3 / 8 and j c oc t 1 / 8 . As a consequence, the ac- 
celerated electron minimum random Lorentz factor and the turbulent magnetic 
field also decrease as inverse power-laws in time (Meszaros and Rees 1993, Rees 
and Meszaros 1992). This implies that the spectrum softens in time, as the syn- 
chrotron peak corresponding to the minimum Lorentz factor and field decreases, 
leading to the possibility of late long wavelength emission (Sari et al. 1998). 

The relativistic expansion is then gradually slowed down, and the blast wave 
evolves in a self-similar manner with a power-law lightcurve. This phase ends 
when so much mass shares the energy that the Lorentz factor drops to 1 (Ayal 
and Piran 2001, Ramirez-Ruiz and MacFadyen 2008). Obviously, this happens 
when a mass E /c 2 has been swept up. This sets a non-relativistic mass scale: 
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For comparison, the Sedov radius of a supernova that ejects a 10 M 0 envelope 
could reach ~ 5 pc or more. 

In GRB sources, with jets which we believe to be highly relativistic, the orienta- 
tion of the jet axis with respect to our line of sight will strongly affect the source’s 
appearance (Dernier 1995, Dalai et al. 2002, Granot et al. 2002, Ramirez-Ruiz 
et al. 2005b, Granot 2007), because radiation from jet material will be Doppler 
beamed in the direction of motion. Attempts to understand the luminosity func- 
tion of GRBs may have to take into account the statistics of orientation, collima- 
tion,and velocity of the jet, as well as the jet’s intrinsic radiation properties. 
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Although our proposed synthesis of GRB physical properties is highly conjectural 
and far from unique, we hope it will provide a framework for discussing the 
integrated properties of these objects. We conclude the review by discussing 
how future observations, experiments and theoretical studies should enhance our 
understanding of the physical properties underlying GRBs. 

VI Prospects 

Gamma-ray burst studies, especially the afterglow-enabled studies of the last ten 
years, remain a young field. The years ahead are thus sure to bring astonishing 
discoveries as the capabilities and experience of observers improve, theorists make 
more and stronger ties to physical theory, and new and upgraded facilities open 
vistas. In this section, we summarize the instrumental capabilities and theoretical 
opportunities for near-term progress in GRB research. 

VI. 1 Facilities 

Table 2 summarizes recent and near-term future GRB missions, planned and 
proposed. Swift (Gehrels et al. 2004) is now the primary mission and has excellent 
prospects for continued operation, with an orbit that will be stable until at least 
2020. Fermi and AGILE, also in continuing operations, are providing added burst 
detections with simultaneous high energy (>100 MeV) coverage that promises to 
redefine the maximum energies and Lorentz factors that GRB engines are capable 
of producing. 

The SVOM mission (Paul et al. 2008), currently under development, promises 4 to 
300 keV coverage and Swift - like slews that will bring X-ray and optical telescopes 



f)§ bear on burst positions. The proposed JANUS small explorer (Rorning 2008) 
would focus on high-redshift bursts, detecting prompt emission over the 1 to 
20keV band and slewing to observe afterglows with a near-infrared telescope 
(50cm aperture, 0.7 to 1.7 /jrn coverage). The proposed large mission EXIST has 
also added fast-response slews and focusing X-ray and near-infrared telescopes to 
its original complement of hard X-ray imaging detectors (Grindlay 2008). 
Current and near-future ground-based facilities include a fast-growing array of 
robotic telescopes primed to respond to burst alerts, ongoing improvements to the 
instrumentation and capabilities of large-aperture telescopes, a new generation 
of air-Cerenkov TeV facilities, and with the EVLA initiative, the most significant 
upgrade to a high-sensitivity radio facility in decades. 

Rapid follow-up of GRB discoveries is occurring both in space and on the ground; 
the Swift UVOT observes ~80 GRBs per year within 2 minutes, detecting about 
half of them. On the ground, fast new telescopes such as ROTSE-III (Akerlof 
et al. 2003), RAPTOR (Wozniak et al. 2006), Pi-of-the-sky (Burd et al. 2005), 
and REM (Zerbi et al. 2001) are able to point at GRBs within 10 — 20 s, and 
larger facilities like the seven-band GROND (Greiner et al. 2008) respond on 10- 
minute timescales. Ready availability of sensitive CCDs and HgCdTe detectors, 
combined with rapid-slew mounts and autonomous software systems, should en- 
able further expansion of rapid- response telescopes with more larger ( D ^ 2 m) 
facilities anticipated in the near future. 

GRB-related programs continue to compete successfully for time on premier op- 
tical facilities. Large-aperture ( D ;> 6 m) telescopes provide the spectroscopic 
observations necessary for GRB redshift measurements. At the VLT, a rapid 
response mode has gone further to provide time-sequence observations revealing 



variable absorption from the host galaxy, and the impending commissioning of 
the X-Shooter spectrograph (Kaper et al. 2009) will enable full UV to NIR charac- 
terization of afterglow spectra with a single integration. In addition to the TOO 
opportunities provided by GRB alerts, multiple host galaxy survey programs are 
underway. 

At radio wavelengths, the workhorse facilities have been the VLA (e.g., Chandra 
et al. 2008) and WSRT (e.g., van der Horst et al. 2007), providing data primarily 
in the 1 — 10 GHz at sensitivities in the 0.2 — 1.0 nrJy range; the Expanded VLA 
(EVLA) upgrade will improve the sensitivity of that facility to ~10 /iJy. Looking 
ahead, by 2012, the ALMA array will be on line operating in the higher frequency 
range 90 — 950 GHz with > 100 times the sensitivity of the VLA. The peak in 
the synchrotron spectrum for a wide range of GRB parameters lies in the ALMA 
range, and should make it a powerful future tool for radio observations. 
Observations of GRBs at TeV energies can be performed by both narrow-field air 
Cerenkov facilities and wide- field water Cerenkov detectors; the latter approach, 
at Milagro, yielded the tentative detection of prompt TeV emission for GRB 
970417A (Atkins et al. 2003). There has been no detection to-date of prompt or 
afterglow emission with the air Cerenkov facilities, but significant effort in pur- 
suit of burst alerts is underway at MAGIC, HESS, and VERITAS. The MAGIC 
dish, in particular, has a rapid-response mode that has provided observations of 
multiple bursts within a minute of trigger (Albert et al. 2006). Future facilties 
will seek to lower detection thresholds to E ~ 100 GeV, which would provide a 
significantly-expanded horizon within which GRB sources will be visible, rather 
than attenuated by photon-photon interactions. 

In space, the upcoming HST SM4 promises to revive that facility, not only pro- 



Qftling new and resuscitated instruments, but also the gyros necessary for flexible 
and fast-response scheduling. In the future, the capabilities of JWST will provide 
an excellent resource for high signal-to-noise observations of highly obscured and 
high-redshift GRBs. 

Within the X-ray band, the Swift XRT has redefined all expectations for the char- 
acterization of X-ray afterglows, and at the same time, proved a highly effective 
facility for refining the multi-arcmin localizations provided by other GRB mis- 
sions. Chandra , XMM, and Suzaku continue to be active in GRB observations, 
and future missions including the Indian ASTROSAT and Japan/US Astro-H 
promise added capabilities in the near future. Dramatic improvements will await 
the next-generation GRB facilities (Table 2) or the arrival of the International 
X-ray Observatory. 

VI. 2 Multimessenger Aspects 

Given the rapid ongoing expansion of the capabilities of non-electromagnetic de- 
tector facilities, and the extreme luminosity and time-specificity of GRB sources, 
it will not be long before we either have the first coincident, multi-messenger de- 
tection of a GRB, or realize limits on the non-electromagnetic emissions of GRBs 
that challenge our current understanding. 

The same shocks which are thought to accelerate electrons responsible for the 
non-thermal 7 -rays in GRBs should also accelerate protons, leading ultimately to 
copious emission of high-energy neutrinos (Waxrnan 2004b, Waxrnan 2004c, Wax- 
man 2006, Dernier and Holmes 2005). The maximum proton energies achievable 
in GRB shocks are E p ~ 10 20 eV, comparable to the highest energies measured 
with large cosmic ray ground arrays (Hayashida et al. 1999). For this, the accel- 



eration time must be shorter than both the radiation or adiabatic loss time and 
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the escape time from the acceleration region (Waxrnan 1995). The accelerated 
protons can interact with the fireball photons, leading to charged pions, muons 
and neutrinos. For internal shocks producing observed 1 MeV photons this im- 
plies > 10 16 eV protons, and neutrinos with 5% of that energy, e v > 10 14 eV 
(Waxrnan and Bahcall 1997). Another copious source of target photons in the 
UV is the afterglow reverse shock, for which the resonance condition requires 
higher energy protons leading to neutrinos of 10 17 — 10 19 eV (Waxrnan and Bah- 
call 1999). Whereas photon-pion interactions lead to higher energy neutrinos 
and provide a direct probe of the shock proton acceleration as well as of the pho- 
ton density, inelastic proton-neutron collisions may occur even in the absence of 
shocks, leading to charged pions and neutrinos (Derishev et al. 1999) with lower 
energies than those from photon-pion interactions. The typical neutrino ener- 
gies are in the ~ 1-10 GeV range, which could be detectable in coincidence with 
observed GRBs. This is the province of projects like AMANDA, IceCube and 
ANTARES. Neutrino astronomy has the advantage that we can see the universe 
up to ~ EeV energies. By contrast, the universe becomes opaque to 7 -rays above 
~ TeV energies through absorption by the infrared background. 

The last and most challenging frontier is that of gravitational radiation, which is 
largely unknown territory. A time-integrated luminosity of the order of a solar 
rest mass (~ 10 54 erg) is predicted from progenitor models involving merging 
compact objects, while that from collapsar models is less certain, and expected 
to be lower by at least one order of magnitude. 

Ground-based facilities, like LIGO, TAMA and VIRGO, are currently seeking 
the first detection of these stellar-scale, high-frequency {v 50 Hz) sources. The 



iSservation of associated gravitational waves would be facilitated if the merg- 
ers involve observed short GRB sources; and conversely, it may be possible to 
strengthen the case for (or against) NS-NS or NS-BH progenitors of short bursts 
if gravitational waves were detected (or not) in coincidence with some events. 
The technical challenge of achieving the sensitivities necessary to measure waves 
from assured sources should not be understated; neither, however, should the 
potential rewards. 

The Enhanced LIGO interferometers will be on line in 2009 with ability to detect 
NS binary mergers to 20Mpc. The Advanced LIGO interferometers on line in 
2014 will extend the distance to 200 Mpc. Short GRBs, if produced by merg- 
ers as is thought to be the case, predict a cosmological rate density of > 10 
Gpc^ 1 yr -1 with likely rate of 300 Gpc 1 yr 1 (O’Shaughnessy et al. 2005, 
O’Shaughnessy et al. 2007, O’Shaughnessy et al. 2008b, O’Shaughnessy et al. 
2008a, Nakar 2007b). This translates into an Advanced LIGO detection rate 
~10 yr -1 . The density estimates include mergers with 7— ray jets not aimed in 
our direction, so not all LIGO detection would be in coincidence with GRBs. 
LIGO is already providing useful upper limits, as with GRB 070201 described in 
Section III. 

VI. 3 Cosmology 

One of the frontiers of modern cosmology lies at high redshift, z ;> 6, when the 
first non-linearities developed into gravitationally-bound systems, whose internal 
evolution gives rise to stars, galaxies, and quasars; and when the light emitted 
from these first collapsed structures diffuse outwards to reionize the Universe. 
As the (temporarily) brightest source of photons in the cosmos, the demise of 



these first generations of massive stars in GRB explosions defines the challenge of 
elucidating the end of the cosmic “dark ages” (Lamb and Reichart 2000, Bromm 
and Loeb 2002). 

Apart from revealing a site of high-redshift star-formation, each such high-redshift 
burst has the potential to help constrain local element abundances in its host 
galaxy, information that will be impossible to gather by other means until the 
advent of D ^ 20 m telescopes. Even more exciting, each burst has the potential 
to reveal the extent of intergalactic reionization, at its redshift and along that 
line of sight (Miralda-Escude 1998, McQuinn et al. 2008, Mesinger and Furlanetto 
2008). 

At the highest redshifts, z ^ 10, there is growing theoretical evidence indicating 
that the first luminous objects to form were very massive stars M > 100 Mq. 
Depending on whether these stars retain their high masses until death, and 
whether a fast-rotating core is prerequisite to the GRB phenomenon, these Pop III 
stars might provide the progenitors for the most luminous, highest-redshift GRBs 
(Heger et al. 2003). 

The first GRBs and supernovae may also be important for another reason: they 
may generate the first cosmic magnetic fields. Mass loss (e.g. via winds) would 
disperse magnetic flux along with the heavy elements. The ubiquity of heavy 
elements in the Lyman alpha forest indicates that there has been widespread 
diffusion from the sites of these early supernovae, and the magnetic flux could 
have diffused in the same way. This flux, stretched and sheared by bulk motions, 
could be the “seed” for the later amplification processes that generate the larger- 
scale fields pervading disc galaxies. 
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What can we expect in the wave of matching theoretical progress? This is more 
difficult to disucss because theory often develops on a shorter timescale than ob- 
servations and experiments, and so we cannot forsee most future developments. 
Although some of the features now observed in GRB sources (especially after- 
glows) were anticipated by theoretical discussions, the recent burst of observa- 
tional discovery has left theory lagging behind. There are, however, some topics 
on which we do believe that there will be steady work of direct relevance to 
interpreting observations. 

One of the most important is the development and use of hydrodynamical codes 
for numerical simulation of GRB sources with detailed physics input. Exist- 
ing two and three dimensional codes have already uncovered some gas-dynanrical 
properties of relativistic flows unanticipated by analytical models (McKinney and 
Blandford 2008) but there are some key questions that they cannot yet address. 
In particular, higher resolution is needed because even a tiny mass fraction of 
baryons loading down the outflow severely limits the maximum attainable Lorentz 
factor. What is more, jets are undoutedly susceptible to hydrostatic and hydro- 
magnetic turbulence. We must wait for useful and affordable three dimensional 
simulations before we can understand the nonlinear development of instabilities. 
Well-resolved three dimensional simulations are becoming increasingly common 
and they rarely fail to surprise us. The symmetry- breaking involved in tran- 
sitioning from two to three dimensions is crucial and can lead to qualitatively 
new phenomena. A particularly important aspect of this would be to link in a 
self-consistent manner the flow within the accretion disk to that in the outflow- 


ing gas, allowing for feedback between the two components. A self-consistent 



model incorporating inflow and outflow must also explain how some fraction of 
the material can acquire more than its share of energy (i.e. , a high entalphy or 
p/p)- 

Particle acceleration and cooling is another problem that seems ripe for a more 
sophisticated treatment. Few would dispute the statement that the photons which 
bring us all our information about the nature of GRBs are the result of particle 
acceleration in relativistic shocks (Medvedev and Spitkovsky 2008, Spitkovsky 
2008, Ranrirez-Ruiz et al. 2007, Katz et al. 2007) or turbulence (Goodman and 
MacFadyen 2007, Zhang et al. 2009, Couch et al. 2008). Since charged particles 
radiate only when accelerated, one must attempt to deduce from the spectrum 
how the particles are being accelerated, why they are being accelerated, and to 
identify the macroscopic source driving the nricrophysical acceleration process. 
Collisionless shocks are among the main agents for accelerating ions as well as 
electrons to high energies whenever sufficient time is available (e.g., Blandford 
and Eichler 1987, Achterberg et al. 2001). Particles reflected from the shock and 
from scattering centres behind it in the turbulent compressed region have a good 
chance of experiencing multiple scattering and acceleration by first-order Fermi 
acceleration when coming back across the shock into the turbulent upstream re- 
gion. Second-order or stochastic Fermi acceleration in the broadband turbulence 
downstream of collisionless shocks will also contribute to acceleration. In addi- 
tion, ions may be trapped at perpendicular shocks. The trapping is a consequence 
of the shock and the Lorentz force exerted on the particle by the magnetic and 
electric fields in the upstream region. With each reflection at the shock the par- 
ticles gyrate parallel to the motional electric field, picking up energy and surfing 
along the shock surface. All these mechanisms are still under investigation, but 



fSere is evidence that shocks play a most important role in the acceleration of 
cosmic rays and other particles to very high energies. 

Another topic on which further work seems practicable concerns the kinematics of 
ultrarelativistic jets (Granot 2007). While it seems probable that we are using the 
correct ingredients of special relativity and a collimated outflow, it is equally true 
that no detailed model yet commands majority of support. We can still expect 
some surprises from studies related to the appearence of relativistic shocks in 
unsteady jets. 

The most interesting problem remains, however, in the nature of the central 
engine and the means of extracting power in a useful collimated form. In all 
observed cases of relativistic jets, the central object is compact, either a neutron 
star or black hole, and is accreting matter and angular momentum. In addition, in 
most systems there is direct or indirect evidence that magnetic fields are present 
- detected in the synchrotron radiation in galactic and extragalactic radio sources 
or inferred in collapsing supernova cores from the association of remnants with 
radio pulsars. This combination of magnetic field and rotation may be very 
relevant to the production of relativistic jets. 

VII Conclusions 

Thanks primarily to the burst discoveries and observations of the Swift satellite, 
the last five years have been tremendously productive ones for GRB research. The 
identification of short burst afterglows has confirmed long-held suspicions that 
GRBs have at least two fundamentally different types of progenitors; subsequent 
studies of short burst afterglows and host galaxies have furnished valuable infor- 
mation on the nature of their progenitors and provided hopeful indications for 



next-generation gravity- wave observatories. The discovery of the first two bursts 
at z > 6, before sources reionized most of the hydrogen in the universe, has proven 
the value of GRBs as probes of the earliest cosmic epochs and and extended GRB 
observations beyond the redshifts of the most distant known quasars. A flood of 
prompt burst alerts has fed the queues of more than a dozen dedicated robotic 
telescopes, and prompted fast-response multi-epoch high-resolution spectroscopy 
from the largest telescopes. Overall, bright GRBs continue to attract the atten- 
tion of astronomers of all types, with premier facilities across the electromagnetic 
spectrum poised to respond to the next spectacular event. In the near future, we 
hope that gravity-wave and high-energy neutrino astronomers will be rewarded 
for their decades of persistent effort. 

The complexity - not to mention sheer volume - of data in the Swift era have 
inevitably raised challenges to the prior interpretation of bursts and their after- 
glows. However, careful consideration of the biases inherent to the Swift mission 
observing strategies, along with the large number of events now available for 
analysis, are gradually enabling the construction of a single coherent picture, via 
a multidisciplinary approach that addresses data from across the electromagnetic 
spectrum. It is one of the challenges of contemporary research to infer the un- 
derlying physical structure of GRBs in the belief that this is simple, despite the 
complex character of the observations. 

GRBs provide us with an exciting opportunity to study new regimes of physics. 
As we have described, our rationalization of the principal physical considera- 
tions combines some generally accepted features with some more speculative and 
controversial ingredients. When confronted with observations, it seems to accom- 
modate their gross features but fails to provide us with a fully predictive theory. 



What is more valuable, though considerably harder to achieve, is to refine models 
like the ones advocated here to the point of making quantitative predictions, and 
to assemble, assess and interpret observations so as to constrain and refute these 
theories. What we can hope of our present understanding is that it will assist us 
in this endeavour. 

There are good prospects for a continued high rate of discovery going forward. 
Swift is likely to remain operational for many years and Fermi and A GILE, now 
on orbit, are providing new insights at higher energies. Instruments currently 
in development or planning have the potential to enable further qualitative ad- 
vances. GRBs are among the most extraordinary of astronomical phenomena 
and will, with our present and future capabilities, continue providing a unique 
window into the extreme reaches of the Universe. 
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Table 1 Typical Parameters of the Canonical Swift X-ray Lightcurve 


Phase 

Start T 

Decay Index 4- 

Approximate 


(s) 


Frequency 

Steep decline 

10 1 - 10 2 

>3 

50% 

Shallow slope 

10 2 - 10 3 

0.5 

60% 

Classical afterglow 

10 3 - 10 4 

1.3 

80% 

Jet Break late phase 

10 5 - 10 6 

2.3 

20%* 

X-ray flares 

10 2 - 10 4 


50% 


+ decay index a defined by F = F a t~ a 

* Of the 80% with no observed jet break, about half had afterglow observations terminate 
before expected time of jet break 
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Table 2 Recent and Future GRB Missions 
Mission Trigger Energy FOV 

Range 

BATSE 20 keV-1.9 MeV 4 t r sr 

(LAD) 

10 keV-100 MeV 


(SD) 


HETE- 

6 - 400 keV 

3 sr 


2 




Swift 

15 - 150 keV 

1.4 sr 


AGILE 

30 MeV - 50 GeV 

~3 sr 


Fermi 

20 MeV - 300 GeV 

> 2 

sr 


(LAT) 

(LAT) 



8 keV to 1 MeV 

9.5 

sr 


(GBM - LED) 

(GBM) 



150 keV to 30 MeV 




(GBM - HED) 



SVOM 

4 keV - 300 keV 

2 sr (CXG) 


(CXG) 




50 keV-5 MeV 

X 

o 

00 

89° 


(GRM) 

(GRM) 


JANUS 

1-20 keV 

4 sr 


EXIST 

5 - 600 keV (HET) 

~ 3.6 

sr 


(HET) 


Detector Area Other Wavelengths GRB rate (yr 1 ) 


2025 cm 2 per 300 

LAD 

127 cm 2 per SD 

120 cm 2 X-ray 

5200 cm 2 UV, Opt, X-ray 100 (~ 10% 

SGRBs) 

Hard X-ray 

> 8000 cm 2 50 

(LAT) 

126cm 2 (GBM - 
LED) 

126 cm 2 (GBM - 
HED) 

Opt, X-ray 80 


Near IR, X-ray (high z) 

5.96 m 2 (HET) Opt, near-IR, X-ray 300 
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Figure 1 7 rays are excellent probes of the most energetic phenomena in nature, that 
typically involve dynamical non-thermal processes and include interactions of high energy 
electrons with matter, photons and magnetic fields; high energy nuclear interactions; 
matter- antimatter annihilation and possibly other fundamental particle interactions. 
Here shown are representative spectra vF v oc v 2 N{y) of GRBs (Kaneko et al. 2007; 
Kaneko et al. 2008) together with the Crab pulsar nebula (Kuiper et al. 2001) and the 
galactic black hole candidate Cygnus X-l (McConnell et al. 2002). 
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Figure 2 Apparent isotropic 7 -ray energy as a function of redshift and observed duration. 
The energy is calculated assuming isotropic emission in a common comoving bandpass 
for a sample of short and long GRBs with measured redshifts. This spread in the inferred 
luminosities obtained under the assumption of isotropic emission may be reduced if most 
GRB outflows are jet-like. A beamed jet would alleviate the energy requirements, and 
some observational evidence does suggest the presence of a jet. 
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Figure 3 The evolving synchrotron afterglow of a 7 — ray burst. Shown is a theoretical 
model (Gou et al. 2007) for the afterglow of the Swift GRB 050904. The model is pre- 
sented without extinction and as it would have been observed at redshift z = 2; the burst 
itself occurred at 2 = 6.29. The evolution of the synchrotron peak to lower frequencies 
is clearly visible. More subtle effects including evolution of the synchrotron cooling and 
self-absorption frequencies, and the associated synchrotron self-Compton emission of the 
blastwave at higher frequencies are not readily visible in this model. 
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Figure 4 Duration and redshift distribution for Swift GRBs. Left panel shows the 
duration distribution. The blue histogram is the measured Tqq distribution while the 
orange one is corrected to the source frame: Tgo/(l + z). Right panel shows the redshift 
distribution for Swift GRBs in blue and pre-Swift GRBs in gray. Swift is detecting higher 
redshift bursts on average than pre-Swift. The solid (broad) theory curve illustrates the 
evolution of a comoving volume element of the universe while the dotted (narrow) curve 
is a convolution of the comoving volume with a model for the star formation rate as 
calculated by Porciani and Madau (2001). 
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Figure 5 Schematic diagrams illustrating the most widely discussed correlations between 
various prompt emission properties for long (L), short (S), and underluminous (UL) 
GRBs. (a) Variability scaled to the burst frame versus Ei so (Fenimore and Ramirez-Ruiz 
2000, Reichart et al. 2001, Schaefer 2006). The variability is a measure of the spikiness 
of the lightcurve and is defined as the mean square of the time signal after removing 
low frequencies by smoothing, (b) Spectral lag scaled to the burst frame versus peak 
luminosity (Norris and Bonnell 2006; Gehrels et al. 2006). (c) Speak scaled to the burst 
frame versus Sj so (Amati et al. 2002 for BeppoSAX GRBs; Amati 2006 for Swift GRBs; 
Lloyd-Ronning and Ramirez-Ruiz 2002 for BATSE events), (d) Duration scaled to the 
burst frame versus Ei so . (e) Spectral hardness versus observed duration (Kouveliotou 
et al. 1993). (f) Pulse rise time versus its decay time (Norris et al. 1996). 
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Figure 6 Representative examples of X-ray afterglows of long and short Swift events with 
steep-to-shallow transitions (GRB 050315, 050724), large X-ray flares (GRB 050502B, 
050724), fastly declining (GRB 051210) and gradually declining (GRB 051221a, GRB 
050826; flux scale divided by 100 for clarity) afterglows. 
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Figure 7 X-ray and optical lightcurves of GRB afterglows in the Swift era. Top panel: 
X-ray lightcurves of Swift burst afterglows. Data for long-duration bursts with known 
redshifts, from GRB 050126 to GRB 070724A, have been gathered from the Swift XRT 
lightcurve and spectral data depository at the UK Swift Science Data Centre (Evans et al. 
2007; Evans et al. 2008). They are corrected uniformly to unabsorbed luminosity over 
1.0 — 30.0 keV in the burst rest-frame, using the time-average afterglow spectrum, and 
plotted as a function of rest-frame time (pink lines) . Separately, afterglow lightcurve fits 
(Racusin et al. 2008b), which exclude flaring intervals, are used to construct minimum 
and maximum envelopes (black lines) and confidence intervals (gray bands) on the X-ray 
luminosities of the bursts as a function of rest-frame time: light gray regions delimit bands 
of 10% to 90% confidence, dark gray regions delimit bands of 25% to 75% confidence, 
and the median burst luminosity at any given time is shown (middle black line). Bottom 
panel: optical lightcurves of Swift burst afterglows. Data for long-duration bursts with 
known redshifts and at least “bronze” quality published optical data (Kann et al. 2007), 
from GRB 050408 to GRB070612A. They are corrected uniformly to rest-frame CZ-band 
luminosity using the inferred i?-band (z = 1) light-curves from Kann et al. (2007), 
and plotted as a function of rest-frame time (pink lines). Interpolated and “best fit” 
extrapolated lightcurves are used to generate minimum and maximum envelopes and 
median luminosity estimates (black lines) and confidence intervals (gray regions), as in 
the previous panel. 
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Figure 8 The growth of GRB 030329 with time as measured using VLBI (by Pihlstrom 
et al. 2007 and references therein). In the background are the images from 2003 April 
(15 GHz), 2003 June (8 GHz), 2003 November (8 GHz), and 2005 June (5 GHz) with 
the same intensity scale. The resolution for the four images is not constant in time, but 
is accounted for in the analysis of the source size. 
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Figure 9 Isotropic-equivalent luminosity of GRB X-ray afterglows scaled to t = 11 hr 
(5 keV source frame) after the burst trigger as a function of their isotropic 7 -ray energy 
release (adapted from Nysewander et al. 2008). 
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Figure 10 The radiated energy inventory of Swift GRBs. Left Panel: Summary of the 
isotropic-equivalent total emitted energy of the prompt and afterglow emission for Swift 
GRBs in the source frame (adapted from Kaneko et al. 2007). Right Panel Bottom : Dis- 
tributions of cumulative isotropic-equivalent total emitted energy for Swift GRBs. Right 
Panel Top: Comparisons of collimation-correctecl total emitted y-ray energy, .Ej so (fl/47r), 
of pre-Swift GRBs, where Ei so is used as an upper limit for GRBs with no jet angle con- 
straints. 
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Figure 11 A selection of the host galaxies of long-duration (top row) and short-duration 
(bottom row) gamma-ray bursts, as imaged by the Hubble Space Telescope. An attempt 
has been made to choose pairs of long and short burst host galaxies with comparable 
redshifts; lower-redshift hosts are emphasized as these reveal their structure more readily 
in short exposures. Images are oriented with North up and East to the left, and the 
physical length scale for a one-arcsecond angular distance is indicated in each panel 
(except for GRB 051227); arrows point to the location of the burst where this is known 
to ~pixel precision. Individual burst notes: GRB 030329 was the first classical long GRB 
to be associated with a well-observed spectroscopic supernova (Stanek et al. 2003, Hjorth 
et al. 2003); XRF 020903 was the first X-ray Flash event to yield a redshift measurement 
Soderberg et al. (2004); GRB 050709 was the first short burst with optical afterglow - 
indicated by the arrow detected Fox et al. (2005), Hjorth et al. (2005b); GRB050509B 
was the first short burst with detected afterglow (Gehrels et al. 2005; Bloom et al. 2006) 
GRB 051227 has a faint candidate host, of unknown redshift probably greater than 1, 
visible at the optical afterglow location; the spiral galaxy to the east has redshift z = 0.714 
Foley et al. (2005). Long-burst host images from Wainwright et al. (2007); short-burst 
host images from Fox et al. (2005) and this work. 
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Figure 12 Keck/LRIS spectrum of the afterglow for GRB 080607 (Prochaskaet al. 2009). 
The red dashed lines indicate a model of the intrinsic afterglow spectrum reddened heavily 
by dust in the host galaxy (rest- frame Ay « 3.2). At A « 4900 A one identifies a damped 
Lya profile associated with H I gas near the GRB. The shaded region overplotted on the 
data corresponds to an H I column density Nh = 10 22 ' 7±o l5 cm— 2. The model (green 
solid line) includes absorption from H 2 Lyman-Werner transitions. The line opacity at 
A > 5500 A is dominated by metal-line transitions from gas in the host galaxy and 
includes bandheads of the CO molecule. Suprisingly, this is the only sightline to date 
to show strong molecular absorption (Tumlinson et al. 2007). It also exhibits a roughly 
solar metallicity. The figure also shows the spectral region or features corresponding to 
intergalactic Lya and Mgll absorption. The redshift for GRB 080607 is z = 3.036. 
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Figure 13 Metallicity as a function of 5-band absolute magnitude for the host galaxies 
of short (red) and long (blue) GRBs. The yellow bars mark the 14 — 86 percentile range 
for galaxies at z ~ 0.1 from the Sloan Digital Sky Survey (Tremonti et al. 2004), while 
crosses designate field galaxies at z ~ 0.3 — 1 (Kobulnicky and Kewley 2004). Both 
field samples exhibit a clear luminosity-metallicity relation. The long GRB hosts tend 
to exhibit lower than expected metallicities (Stanek et al. 2006), while the hosts of short 
GRBs have higher metallicities by about 0.6 dex and are moreover in excellent agreement 
with the luminosity-metallicity relation. From Berger (2009). 
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Figure 14 Diagram exhibiting GRB activity over successive decades in radius ranging 
from 10 6 cm to lpc. The first task in attempting to construct a general scheme of GRBs 
is to decide which parameters exert a controlling influence upon their properties. Mem: 
The GRB nucleus (BH or NS) and its magnetosphere. lOOMcm: The accretion flow 
is likely to be embedded in a very active corona. We expect coronal arches, as well 
as large magnetic structures, to be quite common and to be regenerated on an orbital 
timescale. Relativistic outflow from the black hole is proposed to be focused into two 
jets. lOGcm: Even if the outflow is not narrowly collimated, some beaming is expected 
because energy would be channelled preferentially along the rotation axis. The majority 
of stellar progenitors, with the exception of some very compact stars, will not collapse 
entirely during the typical duration of a GRB. A stellar envelope will thus remain to 
impede the advance of the jet. Tcm: This is the typical size of an evolved massive 
star progenitor. A thermal break-out signal should precede the canonical, softer 7 -rays 
observed in GRBs. lOOTcm: Velocity differences across the jet profile provide a source of 
free energy. The most favourable region for shocks producing highly variable 7 -ray light 
curves is above the baryonic or pair-dominated photosphere. lOPcm: The external shock 
becomes important when the inertia of the swept-up external matter starts to produce 
an appreciable slowing down of the ejecta. lEcm: Finally, we come to the end of the 
relativistic phase. This happens when the mass Ej(? has been swept-up. 
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Figure 15 A 360° vista showing the entire sky, with visible structures stretching back in 
distance, time and redshift. The most distant light we observe comes from the radiation 
left over from the Big Bang: the CMB. As we descend the chart, we find the most 
distant objects known, followed by a web of Sloan Digital Sky Survey (SDSS) quasars 
and galaxies. Closer to home, we start to see a collection of familiar “near” galaxies 
(triangles). Also marked are all Swift GRBs with known distances (blue stars); SN 
1997ff, the most distant type la supernova a,t z = 1.7; and the archetypal large galaxy 
cluster, the Coma cluster. The redshift distances of most distant GRBs are comparable 
to the most distant galaxies and quasars (adapted from Ramirez-Ruiz (2006a).) 
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Figure 16 High-z record holders. The history of most distant objects detected in cate- 
gories of galaxies, quasars, and GRBs. (N. Tanvir, priv. comm.) 
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